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ABSTRACT 


An  investigation  of  the  flow  field  generated  by  a 
rearward  facing  step  at  supersonic  speeds  has  been  con¬ 
ducted.  Measurements  were  made  and  evaluated  in  the  on¬ 
coming  boundary  layer,  in  the  free  shear  layer,  the  lip 
shock  and  in  the  wake  (or  reattachment)  shock  region. 

Pitot  and  static  pressure  surveys  were  conducted,  along 
with  hot  wire  measurements  and  flow  visualization  studies. 
The  measured  free  shear  layer  profiles  were  compared  at 
different  free  stream  Mach  numbers,  Reynolds  numbers  and 
step  heights. 

A  new  technique  was  established  for  reducing  the  hot 
wire  data.  This  was  based  upon  the  heat  transfer  correla¬ 
tion  formula,  and  allows  a  broader  range  of  applicability; 
in  both  the  low  Reynolds  number  region  and  in  the  free 
shear  layer  where  the  flow  field  is  changing  rapidly.  This 
technique  when  applied  to  the  measurements  allowed  the 
flow  structure  in  the  recirculating  region  to  be  explored. 

The  low  base  pressure  generated  by  the  step  influences 
the  flow  field  ahead  of  the  step  and  the  extent  of  this 
upstream  influence  has  been  discussed.  A  new  base  pressure 
correlation  parameter  has  also  been  developed  for  fully 
laminar  and  transitional  flows. 
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CHAPTER  I 


INTRODUCTION 

It  has  been  known  for  sometime  that  the  development 
of  the  entire  wake  flow  generated  by  a  blunt-base  slender 
body  at  supersonic  and  hypersonic  speeds  depends  to  a  great 
extent  on  the  characteristics  of  the  near  wake  [1,  2].^ 
Thus,  the  near  wake  problem  has  been  the  subject  of  inten¬ 
sive  studies  in  recent  years.  The  importance  of  the  de¬ 
tailed  understanding  of  near  wake  for  far  wake  studies  can 
be  illustrated  by  the  fact  that  the  separated  bound¬ 
ary  layer  cannot  be  entrained  by  the  recirculating  region 
and  must  pass  downstream  into  the  far  wake.  A  description 
of  the  flow  region  is  shown  in  Figure  1. 

The  near  wake  is  in  essence  a  strong  interaction 
problem.  For  example,  the  base  pressure  determines  the 
degree  of  expansion  and  separation  for  the  approaching 
boundary  layer,  while  the  separation  process  and  the  free 
shear  layer  growth  provide  the  conditions  before  the  recom¬ 
pression  process,  which,  in  turn,  determines  the  base 
pressure,  the  free  shear  layer  controls  the  structure  of 
the  recirculating  region.  Therefore,  the  free  shear  layer, 


1Numbers  in  brackets  refer  to  similarly  numbered 
references  in  the  bibliography. 
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Figure  1.  Flow  Pattern  for  Supersonic  Flow  over  a  Rearward  Facing  Step. 


the  recirculation  and  the  neck  regions  are  intimately 
coupled  with  one  another.  The  final  solution  to  this  prob¬ 
lem  depends  on  the  balance  of  these  regions. 

The  fundamental  difficulty  in  near  wake  problems 
lies  in  the  fact  that  the  well  established  asymptotic  tech¬ 
niques,  that  form  the  backbone  of  fluid  mechanics,  cannot 
be  applied  with  confidence.  For  instance,  the  boundary 
layer  theory  cannot  be  applied  in  the  expansion  corner  and 
the  recirculating  region.  Only  the  Navier-Stokes  equa¬ 
tions  are  adequate  to  describe  the  complete  flow  field. 
However,  a  successful  numerical  solution  to  these  equa¬ 
tions  seems  not  to  be  feasible  in  the  near  future.  New 
avenues  must  be  attempted  in  establishing  a  physically 
realistic  model.  It  is,  therefore,  very  necessary  to  ex¬ 
perimentally  describe  the  detail  flow  structure  in  order 
to  determine  the  relative  importance  of  these  flow  regions 
and  to  verify  the  existing  theoretical  models. 

It  was  Chapman  [3]  who  first  recognized  the  impor¬ 
tance  of  the  free  shear  layer  shed  from  the  body  surface. 

He  introduced  the  concept  of  "dividing  streamline"  sepa¬ 
rating  the  recirculating  fluid  from  the  fluid  entrained  by 
mixing  with  the  external  uniform  inviscid  flow.  However, 
the  uniform  parallel  inviscid  flow  is  a  severe  idealization 
of  the  corner  expansion  process,  which  is  valid  only  when 
the  oncoming  boundary  layer  thickness  is  very  thin  or  the 
free  stream  Reynolds  number  is  very  large. 
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Crocco  and  Lees  [4]  proposed  the  concept  that  the 
compression  in  the  near  wake  at  supersonic  speed  is  a 
smooth  viscous-inviscid  interaction  phenomenon  rather  than 
a  sudden  isentropic  process  as  assumed  by  Chapman.  The 
mixing  between  the  internal  dissipative  flows  and  the  ex¬ 
ternal  nearly  isentropic  streams  is  regarded  as  a  funda¬ 
mental  physical  process.  Their  theory  shows  the  proper 
variation  of  base  pressure  with  Reynolds  number;  however, 
in  order  to  bring  out  the  main  aspects  of  their  theory,  the 
treatments  of  the  recirculating  flow  and  the  corner  ex¬ 
pansion  process  were  somewhat  crude. 

The  problem  of  upstream  influence  of  the  low  base  pressure 
has  been  investigated  by  Scherberg  &  Smith  [18],  Baum  [5],  Olsson 
and  Messiter  [6],  and  Weiss  and  Nelson  [7],  The  problem  of 
expansion  and  separation  of  the  approaching  boundary  layer 
in  the  neighborhood  of  a  sharp  corner  at  the  base  has  been  studied 
by  Scherberg  &  Smith  [18],  Hama  [8],  Weiss  and  Weinbaum  [9],  Adam¬ 
son  [10],  and  Zakky  [11].  The  recirculating  region  has  been  attacked 
by  Batchelor  [12],  Burggraf  [133,  Weiss  [14],  Scherberg 
[15],  Roache  [16],  and  Erdos  and  Zakky  [17] .  Despite  the 
strong  interacting  (or  coupling)  nature  of  the  near  wake, 
so  far,  the  problems  on  separation  corner,  free  shear  layer 
and  recirculating  region  have  been  mainly  analytically  ex¬ 
amined  independently  of  each  other.  The  application  of 
inviscid  theory  to  determine  the  expansion  around  a  separa¬ 
tion  corner  remains  to  be  justified  by  experiments.  There 
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appears  to  be  a  controversy  [8]  as  to  cause  of  the  lip 
shock.  The  effects  of  the  expanded  boundary  layer  on  the 
development  of  the  inner  free  shear  layer. and  the  structure 
of  the  recirculating  region  have  not  been  investigated 
rigorously.  No  attempt  has  been  made  to  correlate  the 
base  pressure  for  the  transitional  flow  up  to  date.  Many 
valuable  experimental  investigations  in  the  near  wake  of 
a  blunt-based  slender  body  have  been  conducted  by  Scherberg 
and  Smith  [18],  Batt  and  Kubota  [19],  Roshko  and  Thomke 
[20],  Lewis  and  Behrens  [21],  Todisco  and  Pallone  [22], 

Rom  [23],  Charwat  and  Yakura  [24],  and  many  others.  How¬ 
ever,  more  information  is  still  needed  in  the  regions  of 
the  upstream  influence,  corner  expansion,  free  shear  layer 
and  the  recirculating  flow. 

The  objectives  of  this  investigation  are  to  study 
experimentally : 

1.  The  upstream  influence,  the  corner  expansion 
process,  and  the  comparisons  with  the  character¬ 
istics  solution  and  with  the  Prandtl-Meyer 
theory . 

2.  The  basic  characteristics  of  the  recirculating 
region. 

3.  The  effects  of  Reynolds  number,  Mach  number  and 
step  height  on  the  separated  boundary  layer. 

4.  The  geometry  of  the  reattachment,  or  wake,  shock. 

The  upstream  influence  of  the  low  base  pressure  was 

studied  through  the  use  of  very  densely  distributed  surface 


pressure  measurements  and  pitot  surveys  of  the  approaching 
layer.  Very  extensive  pitot  surveys  were  conducted  immedi¬ 
ately  downstream  of  the  corner  to  study  the  rapid  expansion 
process  and  the  formation  of  the  lip  shock.  The  free  shear 
layer  growth  was  also  investigated  by  pitot  pressure  meas¬ 
urements.  Because  much  experimental  information  has  been 
accumulated  at  M  =  3.5,  Re  T  =  1.0  x  106,  the  hot-wire 

co  oOjXj  7 

measurements  were  conducted  at  this  condition. 

Hot-wire  anemometer  probes  were  used  to  obtain  the 
absolute  velocity  and  total  temperature  in  the  recircu¬ 
lating  region  and  the  free  shear  layer.  The  Reynolds  number 
based  on  wire  diameter  is  not  only  very  low  in  the  recir¬ 
culating  region  but  also  changes  very  abruptly  across  the 
free  shear  layer.  Because  of  the  wide  range  of  applica¬ 
tion  (M,  ReQ)  and  the  complicated  heat  transfer  formula, 
it  is  almost  impossible  to  apply  the  conventional  hot¬ 
wire  technique.  Therefore,  a  simplified  data  reduction 
method  was  developed  for  this  investigation. 

This  experimental  investigation  is  part  of  the 
Separated  Flow  Program,  Aeronautical  Research  Laboratory, 
(ARL) ,  Thermodynamics  Division,  Wright-Patterson  Air  Force 
Base,  Ohio,  and  a  continuation  of  the  research  efforts  of 
Rom  [23]  and  Scherberg  and  Smith  [18],  The  test  model, 
pitot-probe  and  hot-wire  holders  were  designed  and  built 
at  the  ARL.  Many  pitot-pressure  surveys  analyzed  in  this 
study  were  also  conducted  by  Smith  [25]  during  the  previous 
test  program.  Some  results  of  this  investigation  were 
presented  by  Wu,  Su  and  Scherberg  [26], 
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CHAPTER  II 


EXPERIMENTAL  TECHNIQUES 
1.  WIND  TUNNEL  DESCRIPTION 

Tests  were  conducted  in  Tunnels  A  and  D  of  the  von 
Karman  Gas  Dynamics  Facility  (VKF)  at  the  Arnold  Engineer¬ 
ing  Development  Center  (AEDC)  [27].  Tunnel  A  is  a  con¬ 
tinuous,  closed-circuit,  variable  density  supersonic  wind 
tunnel  with  an  automatically  driven,  flexible  plate-type 
nozzle  and  a  40  by  40  inch  test  section.  The  tunnel  can 
be  operated  at  Mach  numbers  from  1.5  to  6  and  unit  Reynolds 
numbers  per  inch  between  0.025  x  10 6  and  0.59  *  10®. 

Tunnel  D  is  an  intermittent  supersonic  wind  tunnel 
with  a  manually  adjusted,  flexible-plate-type  nozzle  and  a 
12  by  12  inch  test  section.  The  tunnel  operates  at  Mach 
numbers  from  1.5  to  5.0,  at  stagnation  pressures  from  5  to 
60  psia  and  at  stagnation  temperatures  up  to  about  80°F. 

2.  WIND  TUNNEL  MODELS 

The  model,  for  Tunnel  A  is  shown  in  Figure  2.  It  is 
a  two-dimensional  rearward  facing  step  model,  9  inches  wide 
with  variable  step  heights  of  0.443  and  0.75  inches  and 
it  is  the  same  model  as  used  in  the  earlier  tests  [25,  28,  29, 
30].  Sideplates,  fitted  with  optical  quality  glass  windows, 
were  attached  to  the  model  to  obtain  approximately 
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Pressure  Orifice  Location 


Figure  2.  Model  Details 


two-dimensional  flow.  The  diameter  of  the  leading  edges  of 
the  model  and  sideplates  were  less  than  0.002  inch.  The 
maximum  Reynolds  number  based  on  tunnel  free  stream  con¬ 
ditions  and  the  leading  edge  diameter  was  1000.  Eleven 
pressure  orifices  were  located  on  the  model  top  surface, 

34  on  the  lower  reattachment  wall  and  seven  on  the  verti-  ■ 
cal  step  face,  dimensions  of  pressure  orifices  are  speci¬ 
fied  in  Appendix  G.  A  longitudinal  vertical  splitter  plate 
instrumented  with  22  orifices  was  used  for  static  pressure 
measurements  in  and  near  the  recirculating  region,  the 
details  are  shown  in  Figure  3. 

The  model  for  the  Tunnel  D  testing  was  furnished 
by  VKF ,  AEDC.  This  model  spanned  the  test  section  and 
was  12  inches  wide  and  did  not  have  sideplates.  The 
height  of  the  step  could  be  adjusted  to  0.75  inch  and 
0.225  inch.  A  splitter  plate  was  used  for  oil  flow 
studies  and  static  pressure  measurement  [28], 

3.  PITOT  PRESSURE  MEASUREMENTS 

The  model  boundary  layer,  expansion  region  and 
free  shear  layer  were  surveyed  with  a  flattened  combina¬ 
tion  pitot-thermocouple  probe.  The  probe  is  0.01  inch 
high  and  with  a  0.004  inch  by  0.06  inch  opening,  see 
Figure  4.  The  probe  was  always  pointed  upstream  and 
parallel  to  the  free  stream  direction.  No  attempt  has 
been  made  to  correct  the  errors  due  to  flow  inclination 
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All  dimensions  in  inches 


Pressure 


0.750 


Figure  4.  Combination  Thermocouple  Pitot  Probe  Details. 
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and  low  Reynolds  number.  The  height  of  the  probe  center- 
line  can  be  varied  with  an  accuracy  of  ±0.001  inch. 

The  pitot  pressure  measurements  were  made  with  a 
15-psid  transducer  calibrated  for  ranges  of  15,  5,  1,  and 
0.3  psia.  Whenever  possible,  a  near  vacuum  reference  was 
used;  otherwise,  a  2-psid  reference  was  employed.  The 
transducer  is  considered  to  be  accurate  to  within  ±1.0 
percent  of  the  range  used. 

4.  STATIC  PRESSURE  MEASUREMENTS 

Limited  static  pressure  profiles  were  conducted 
with  a  0.25  inch  diameter  "disk  type"  probe.  A  0.047  inch 
diameter  orifice  was  located  at  the  center  of  the  disk  as 
shown  in  Figure  5.  Because  of  the  relatively  large  ori¬ 
fice  diameter,  the  measurements  are  not  accurate  in  the 
regions  of  abrupt  pressure  changes  such  as  the  expansion 
region.  Static  pressure  measurements  were  made  with  a 
5-psid  transducer  reference  to  a  near  vacuum.  The  trans¬ 
ducer  was  calibrated  for  ranges  of  2.0,  1.0,  and  0.3 
psia,  and  was  considered  to  be  accurate  to  within  ±1.0 
percent  of  the  range  used. 

Model  surface  pressures  were  measured  with  a  pres¬ 
sure  scanning  system  which  used  15-psid  or  5-psid  trans¬ 
ducers  reference  to  a  near  vacuum.  The  15-psid  trans¬ 
ducers  were  calibrated  for  ranges  of  5  and  1  psia,  and 
were  considered  accurate  to  within  ±0.4  and  ±1.0  percent. 
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respectively,  of  the  full  scale  of  the  range  used.  The 
lower  surface  pressures,  such  as  base  pressure,  were  meas¬ 
ured  with  the  pressure  scanning  system  using  5-psid  trans¬ 
ducers  reference  to  a  near  vacuum.  The  transducers  were 
calibrated  for  ranges  of  5,  1.67,  and  0.33  psia.  These 
transducers  were  considered  accurate  to  within  ±0.2,  ±0.3, 
and  ±1.0  percent,  respectively,  of  the  range  used. 

5 .  HOT-WIRE  MEASUREMENTS 

The  hot-wire  holders  as  shown  in  Figure  6  were 
designed  so  that  the  measuring  wire  was  offset  from  the 
supporting  stem  which  was  made  [25]  long  and  thin  to 
cause  minimum  disturbance  to  the  flow  structure.  Each 
probe  consisted  of  platinum  -  10  percent  rhodium  wire 
approximately  0.00022  inch  in  diameter  and  0.065  inch  in 
length.  A  thermocouple  was  attached  to  one  needle  of  each 
probe  approximately  0.1  inch  from  the  tip  to  measure  the 
temperature  of  the  wire  support.  All  wires  were  annealed 
and  calibrated  following  the  method  outlined  by  Dewey  [31]. 
From  these  calibration  measurements,  the  linear  tempera¬ 
ture  resistance  relation  was  obtained  as  shown  in  Figure  7. 

The  hot-wire  measurements  were  made  using  constant 
current  equipment  of  Shapiro-Edwards  design.  Wire  resis¬ 
tance  was  measured  by  a  precision  bridge  with  an  accuracy 
of  ±.01  ohm.  Wire  heating  currents  were  set,  using  a 
precision  potentiometer  with  ranges  of  0.1  amp  and  0.01 
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-Wire  Anemometer  Probe  Details. 
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Calibration 


amp  full  scale  with  an  accuracy  of  ±0.1  percent  of  full 
scale.  The  measurements  of  wire  support  temperature  were 
accurate  to  within  ±1.0  degree. 
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CHAPTER  III 


IMPROVED  HOT-WIRE  TECHNIQUE  FOR  LOW  REYNOLDS 
NUMBER  COMPRESSIBLE  FLOWS 

1.  GENERAL  DISCUSSION  ON  HOT-WIRE  TECHNIQUES 

The  hot-wire  methods  have  been  used  for  a  number  of 
years  to  measure  the  mean  or  fluctuating  flow  properties 
of  a  fluid  stream.  The  basic  principle  of  the  hot-wire 
method  may  be  stated  as  follows.  If  a  small  wire  is  placed 
in  a  moving  medium  and  heated  to  a  temperature  higher  than 
that  of  the  medium,  there  will  be  heat  transfer  from  the 
wire  to  the  flow.  The  rate  of  heat  exchange  depends  on 
the  geometrical  and  physical  properties  of  the  wire,  the 
physical  characteristics  of  the  medium,  the  Joule  heating, 
and  the  flow  conditions.  From  the  measured  rate  of  heat 
exchange,  one  can  determine  the  flow  condition.  The  funda¬ 
mental  information  needed  in  order  to  use  a  hot-wire  probe 
is  the  governing  law  of  heat  transfer  from  a  single  fine 
wire  under  the  conditions  of  which  the  measurements  will 
be  performed. 

The  dependence  of  the  Nusselt  number  on  the  Reynolds 
number  for  incompressible  flow  was  first  established  by 
King  [32].  He  found  that  the  Nusselt  number,  Nu,  is  a 
linear  function  of  the  square  root  of  the  Reynolds  number. 
Re ,  i .  e . , 
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Nu  =  A  +  B  Re 


(1) 


1/2 


where  A  and  B  are  constants.  This  equation  may  be  reduced 
to: 


Im  Rm 
m  m 

R_  -  R 
m  a 


=  K'  +  B"V 


1/2 


(2) 


where  I  is  the  current  passing  through  the  wire  at  the 

wire  resistance  R„,,  R,  is  the  wire  resistance  at  zero 

m  a 

current.  A'  and  B “  are  constants  and  V  is  the  velocity  of 
the  flow. 

From  Equation  (2),  one  knows  that  once  A'  and  B' 

are  experimentally  calibrated,  the  velocity  V  can  be  deter- 
2 

mined  if  (Im  Rm)/ " Ra)  is  measured.  For  compressible 
flow.  Equation  (1)  takes  a  more  complicated  form  [33] , 
i.e. , 

Nu  =  A  +  B(M)  Re1/2  (3) 

where  B,  previously  a  pure  constant,  is  now  a  function  of 
the  flow  Mach  number  M.  Consequently,  the  calibration 
process  is  more  complicated.  Now,  if  the  flow  density  is 
very  low  such  that  in  the  region  of  the  slip  or  the  free 
molecular  flows,  the  functional  relationship  between 
Nusselt  and  Reynolds  numbers  becomes  very  complicated. 
Also,  at  very  low  Reynolds  number,  the  velocity  induced  by 
buoyancy  may  be  of  greater  importance.  The  natural  con¬ 
vection  shall  be  included  in  the  heat  transfer  analysis. 
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Neglecting  natural  convection,  Dewey  [34]  gave  the 
following  functional  relationships: 

Nuq  =  A(Reo)  B(ReQ,  M)  ReoN(ReoJ  (4) 

T  =  T  (M,  Kn)  (5) 

aw  aw  '  ' 

where 

Rec  =  Re0(M,  Tq#  pg) 

TTY  1/2  M 

Kn  =  (— )  - 

2  Re„ 

00 

The  detail  expressions  are  very  much  involved,  and 
will  be  discussed  later.  One  -realizes  that  it  may  not  be 
wise  to  attempt  the  calibration  technique.  A  new  method 
must  be  developed  to  overcome  this  difficulty. 

Most  of  the  hot-wires  are  of  a  rather  fine  wire; 
their  length-to-diameter  ratio  usually  varies  from  100  to 
300.  The  wire  materials,  on  the  other  hand,  are  good  heat 
conducting  metals.  A  substantial  portion  of  the  heat 
generated  in  the  wire  is  conducted  away  toward  the  cooler 
ends  of  the  wire  supports.  In  order  to  obtain  the  correct 
value  of  heat  transferred  by  convection  to  the  surrounding 
air,  this  conduction  heat  loss  must  be  separated  from  the 
total  heat  loss. 

In  the  continuum  flow,  the  conductive  heat  loss  is 
usually  only  a  small  fraction  of  the  convective  heat  loss. 
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if  the  wire  length-to-diameter  ratio  l/d  exceeds  100.  The 
end  loss  can  be  taken  into  account  as  a  correction  factor. 
It  only  modifies  the  total  heat  loss  by  a  constant  factor 
which  is  almost  independent  of  the  flow  velocity.  How¬ 
ever,  in  the  low  density  flows,  the  end  loss  encountered 
may  be  very  appreciable  even  for  a  large  &/d  case.  To 
obtain  a  reliable  and  meaningful  result,  great  care  must 
be  exercised  when  the  end  loss  correction  technique  is 
employed.  A  hot-wire  calibrated  in  a  uniform  flow  has  the 
end  loss  effect  included  in  the  calibration.  However,  if 
the  probe  supports  are  subjected  to  a  large  temperature 
gradient,  the  support  temperature  is  different  from  that 
for  uniform  flow,  and  the  change  in  end  losses  will  in¬ 
validate  a  calibration  made  in  uniform  flow.  It  then 
becomes  necessary  to  calculate  the  end  losses  with  the 
measured  support  temperature. 

Dewey  [31]  developed  an  end  loss  correction  technique 

for  the  hot-wire  measurements  in  low  Reynolds  number 
hypersonic  flow.  The  end  loss  correction  factors,  and 
i|/R  were  derived  in  terms  of  NuQ,  then  ipR  and  the 
empirical  heat  transfer  correlation  formula  were  coupled 
to  solve  for  flow  field  parameters.  Also,  the  expression 
for  ijjN  is  formulated  with  some  approximations.  Because 
the  method  of  successive  substitution  was  used  to  solve 
the  simultaneous  equations,  the  calculation  was  found  to 
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be  nonconvergent  for  M  <  1.5  [35].  Dewey's  technique  is 
reviewed  in  Appendix  A. 

If  one  is  interested  in  the  hot-wire  measurement 
in  the  low  Reynolds  number  supersonic  or  subsonic  flow 
where  the  end  correction  is  even  more  important,  one  must 
keep  in  mind  that  approximations  shall  be  kept  at  a 
minimum  whenever  possible,  and  that  any  approximations 
employed  need  careful  justification. 

One  of  the  objectives  of  this  investigation  is  to 
explore  the  basic  characteristics  of  the  recirculating 
region  by  hot-wire  measurements.  The  recirculating  region 
was  once  considered  a  "dead  air"  region.  The  velocities 
in  a  portion  of  this  region  are  low  subsonic  and  the  Rey¬ 
nolds  number  based  on  the  hot-wire  diameter  is  also  very 
low.  In  order  to  obtain  meaningful  measurements,  a  direct 
and  more  simple  technique  was  developed  to  determine  the 
convective  heat  transfer  coefficient  h  and  the  recovery 
temperature  T  of  the  wire  [36).  This  technique,  of 

aw 

course,  is  applicable  to  the  supersonic  and  hypersonic 
flows  as  well. 


2.  DIRECT  DETERMINATION  OF  h  AND  Taw 

As  discussed  in  the  last  section,  if  one  follows 
Dewey's  approach,  one  must  solve  the  following  two  equa¬ 
tions  simultaneously  for  the  flow  field  parameters: 


22 


t 


r 


.;st-  •. 


Taw(M,Kn)  =.  TaWm  *R  Wu0) 


Nuo(M,Reo)  =  Num  *N  (Kao) 


Because  of  the  complicated  relationship  of  the  above  equa¬ 
tions,  the  method  of  successive  approximation  may  not 
always  produce  convergent  solutions  to  the  system  of  equa¬ 
tions.  Also,  the  approximations  for  i|>N  may  introduce 
error  in  low  Reynolds  number  supersonic  or  subsonic  flow. 

It  will  be  shown  that: 

1.  The  end  loss  correction  factors,  ^  and  \pR,  are 
not  really  necessary  for  the  hot-wire  data 
reduction,  and  the  approximations  required  to 
obtain  can  be  eliminated. 

2.  h  and  T  „  can  be  reduced  directly  from  the  hot- 

aw 

wire  measurements,  i.e., 

h  -  h‘V  V 

Taw  “  Taw^Rm'  *1^ 

3.  The  flow  field  parameters  can  be  solved  much 
easier  from  the  following  two  equations: 

Taw  “  Taw(M'  Kn> 
hd/kD  =  Nu0(M,Re0) 

Parts  1  and  2  above  are  discussed  in  this  section, 
and  part  3  will  be  discussed- in  the  next  section. 

The  steady  state  heat  balance  at  any  section  of  a 
finite  wire  is 
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(6) 


,  ird2£  d2T 

Xl  R  +  (— 1 1  kw  ^  =  *d£h  <T  -  Taw> 

Where  the  right  hand  side  of  the  equation  represents  the 

convective  heat  loss,  the  second  term  on  the  left  hand  side 

2 

is  the  heat  loss  by  conduction,  and  I  R  is  the  Joule  heat- 
ing.  This  equation  involves  the  following  assumptions: 

1.  The  convective  heat  transfer  coefficient,  h, 
is  a  constant  and  is  equal  to  the  value  which 
would  be  measured  by  an  infinite  wire  at  the 
same  current.  This  implies  that  the  external 
flow  is  uniform  along  the  wire. 

2.  The  adiabatic  wire  temperature,  T  ,  is  a  con- 

dw 

stant  and  is  equal  to  the  value  of  an  infinitely 
long  wire. 

3.  Natural  convection  and  radiation  are  negligible. 

4.  The  wire  thermal  conductivity  is  independent 
of  temperature  within  the  heated  range. 

5.  Wire  thermal  resistance  varies  linearly  with 
temperature ,  i . e . , 

R  =  a  +  bT  (7) 

and  R  is  the  resistance  of  the  wire  of  length 

£. 

6.  Each  cross-section  of  the  wire  is  at  uniform  T, 
h  and  Taw. 
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The  heat  transfer  equation  can  be  rearranged  as 


2 

d  T 


dx 


-  m  T  +  n  =  0 


(8) 


where 


n  = 


4(7Td*h  Taw  +  a  I/) 


ir  d*6  i  k 


(9) 


w 


9  4  (ndS-h  -  b  I  ) 

ill 

in  *  . 

it  i  k 

w 


(10) 


and  the  boundary  conditions  are  T  =  Tg  at  x  =  ±1/2  where 
T  is  the  temperature  of  the  supports.  The  solution  to  the 

b 

above  equation  is 


n 


n 


cosh  mx 


T  =  (-j  -  T  >  - 

nr  nr  cosh  mi/2 


(11) 


and  the  measured  resistance  of  the  wire  is 


1/2 


%  - 


R(T>  dx  =  a  +  b  4  -  b("/.m.2.'-,.I-sl.  7  <12> 

nr  mi/2  1 


Note  that  h  and  Taw  are  the  two  unknowns  that  one  wishes  to 

determine  from  the  hot-wire  measurements,  h  and  T  are 

aw 

related  to  the  measurable  R  and  I  by  the  above  nonlinear 

m  m  J 

algebraic  equation.  If  two  wire  resistances  are  measured 
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at  two  different  currents,  two  equations  may  be  obtained 
from  Equation  (12)  to  determine  h  and  T  .  i.e., 

aW 

^  <h*  Taw' 

^2  ~  Rra2  Taw'  ^2^ 

By  operating  at  several  currents,  many  pairs  of 
equations  can  be  solved  for  h  and  T  .  However,  as  ob- 

dw 

served  by  Vrebalovich  [37],  the  solutions  depend  on  the 
degree  of  wire  heating. 

This  problem  can  be  overcome  by  calculating  h  and 
with  the  measurements  corresponding  to  the  no  over¬ 
heating  condition,  i.e.,  I  0.  The  following  two  equa¬ 
tions  can  be  obtained  from  Equation  (12) s 


I  -*■  0 


=  (a  +  b  Taw)  - 


b(Taw  -  Ts)  tanh  S 


(15) 


and  a  plot  of  Rm  versus  I*  allows  the  quantities 


5 

di. 


p. 


m 


I  0 


I  0 


to  be  extrapolated  from  the  measurements  which  are  shown 
in  Figure  8.  Equations  (13)  and  (14)  can  be  solved  for  h 
and  ?aw,  the  details  are  described  in  Appendix  B.  Typical 
results  computed  with  and  without  overheat  are  shown  in 
Figure  9. 

3.  REDUCTION  OF  FLOW  PROPERTIES 


It  is  shown  that  the  hot-wire  is  capable  of  pro¬ 
viding  two  thermodynamic  measurements:  h  and  Taw.  If 
sufficiently  accurate  relations  are  also  known  to  relate 
them  to  the  other  parameters  of  the  flow  (such  as  Mach 
number  and  Reynolds  number) ,  the  flow  field  parameters  can 
be  defined  if  another  thermodynamic  measurement  such  as 
the  static  or  pitot  pressure  is  conducted. 

Because  of  the  wide  range  of  Mach  numbers  and 
Reynolds  numbers  encountered  (from  the  shear  layer  to  the 
recirculating  region) ,  it  is  necessary  to  have  the  ana¬ 
lytical  expressions  for  the  empirical  correlations  suit¬ 
able  for  computer  calculation.  The  correlation  formula 
proposed  by  Dewey  was  employed,  because  the  accuracy  of 
his  expressions  is  superior  to  several  of  the  earlier 
correlation  relations  (34).  The  Nusselt  number  and 
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Figure  8.  Typical  Measurements  in  Recirculating  Region. 


recovery  factor  are  related  to  the  Reynolds  number  and  Mach 
number  by  the  following  expressions: 

Nuo(ReQ/  M>  =  UReo,»)  ip  (ReQ,M)  (16) 

T 

n(M,  Kn)  =  =  n  +  n*(of  -  nj  (17) 

T  c  i  c 

O 

-2  3 

in  the  range  2  x  10  <  ReQ  <  10  and  M  >  0.2,  where 

^(ReQ,a>)  represents  the  heat  transfer  correlation  at  the 
hypersonic  Mach  number  as  a  result  of  the  hypersonic 
"freeze."  n  is  the  ratio  of  wire  recovery  temperature  to 
the  free  stream  stagnation  temperature,  and  the  sub¬ 
scripts  f  and  c  denote  free  molecular  and  high  Reynolds 
number  continuum  limits  respectively. 

The  Nusselt  number  varies  monotonically  from  the 

0  5 

continuum  heat  transfer  law  Nuq  -  ReQ’  at  higher  Reynolds 
numbers  to  Nuq  ~  ReQ  in  free  molecular  flow.  Equation  (16) 
and  (17)  are  shown  in  Figure  1  and  2  of  Reference  [34] . 

The  heat  transfer  correlation  is  accurate  to  ±7  percent, 
and  the  normalized  recovery  factor  n*  is  accurate  to 
±10  percent  for  M  >  0.6  [34] .  Therefore,  the  reduced  flow 
properties  shall  only  be  used  in  a  qualitative  manner. 

The  detailed  expressions  are: 

Re0*7114 

5(Re  .»)  =  Re^fO.  1400  +  0.2302  ( - - - pr-^yr) 
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Knco  =  <A«/d>  *  (if A/2) 


(M/ReJ 


1/2 


Since  p  is  measured,  Equation  (16)  and  Equation  (17)  can 
s 

be  solved  for  M  and  TQ  by  Newton-Raphson  iteration  scheme. 
This  is  discussed  in  Appendix  C. 

All  computations  of  flow  properties  were  made  using 
y  s*  1.4  and  the  power  law  for  viscosity  and  thermal 
conductivity,  or 


T  0.76 
=  (-£) 

T 


:  T  0.86 

°  =  (J9L) 


It  was  found  that  for  M  <  0.3,  there  was  a  numerical  con¬ 
vergence  problem.  This  is  probably  due  to  the  inaccuracy 
of  Dewey's  correlation  formula  in  that  very  low  Mach 
number  range. 
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CHAPTER  IV 

FLOW  UPSTREAM  OF  SEPARATION 

1.  SURFACE  STATIC  PRESSURE  DISTRIBUTION 

Surface  pressure  measurements  upstream  of  the 
separation  corner  are  presented  in  Figures  10,  11,  and 

Vr 

12  for  M^  =  2.5,  3.5,  5.0  and  various  Reynolds  numbers. 

The  pressure  distribution  is  the  combined  results  of  the 
viscous  weak  interaction  and  the  upstream  influence  of 
the  base  pressure.  The  viscous  interaction  effect  is 
pointed  out  by  Baum  [5]  based  on  Hama's  data  [8].  In  order 
to  isolate  the  weak  interaction  effect,  the  surface  pres¬ 
sure  is  'ormalized  by  p^  which  is  the  measured  pressure  ' 
at  the  orifice  furthest  (1.5  inch)  upstream  from  the 
corner. 

A  significant  fraction  of  the  pressure  drop  occurring 
upstream  of  the  corner  can  be  seen.  This  indicates  that 
the  signal  lower  base  pressure  is  transmitted  through 
the  subsonic  portion  of  the  approaching  boundary  layer, 
and  the  expansion  actually  begins  upstream  of  the  corner, 
which,  consequently  introduces  an  appreciably  normal  pressure 
gradient.  At  the  same  free  stream  Mach  number,  the  degree 
of  upstream  influence  gets  larger  as  the  Reynolds  number 
becomes  lower.  The  length  of  upstream  influence  increases 
with  increasing  Mach  number.  These  seem  to  imply  that  the 
♦from  0.2  to  2.0  x  10b 
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Figure  10.  Surface  Pressure  Distributions  at  M( 


boundary  layer  thickness  is  a  characteristic  length  in 
this  region. 


2.  CORRELATION  OF  SURFACE  PRESSURE 

The  surface  pressures  are  correlated  by  the 
parameter  x/6,  where  6  is  the  boundary  layer  thickness  at 
the  separation  corner  computed  as  if  there  were  no  upstream 
influence  [4]. 


6  ,  _ 

-=  [5.0  +  1.2  (Y-l)  M£]//5e^L 
L 


(18) 


The  result  of  correlation  for  free  stream  Mach  number  from 
2.5  to  5.0  and  free  stream  Reynolds  number  per  inch  from 
0.04  x  106  to  0.53  x  io6  is  shown  in  Figure  13.  The 
correlation  reveals  that  the  upstream  influence  is  a 
function  of  6.  The  signal  of  the  lower  base  pressure 
can  be  propagated  upstream  to  a  distance  of  about  five 
times  the  thickness  of  the  approaching  boundary  layer. 

It  appears  from  this  study  that  the  step  height  has  no 
effect  on  the  upstream  influence.  The  results  of  a  6°  wedge 
[8]  are  also  included  in  Figure  13. 

Based  on  Figure  13,  it  is  quite  possible  that, 
right  at  the  corner,  the  pressure  ratio  might  decrease 
to  about  0.528.  This  suggests  that  a  fluid  element 
initially  close  to  the  surface  (just  outside  of  the 
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thin  viscous  layer)  might  accelerate  to  the  sonic  speed 
at  the  corner  as  speculated  by  some  investigators  [6], 

This  will  be  further  discussed  based  on  pitot-probe  meas¬ 
urements  in  the  next  section. 

3.  ORDER  OF  MAGNITUDE  ANALYSIS 

A  typical  pitot  pressure  survey  of  the  approaching 
boundary  layer  at  =  5.0  is  presented  in  Figure  14. 

Because  the  height  of  the  pitot-probe  used  in  this  study 
was  relatively  large  compared  to  the  measured  boundary 
layer  thickness,  there  was  some  probe  interference.  The 
interference  has  been  shown  to  be  severe  at  the  outer  edge 
of  the  boundary  layer  [38,  39,  40].  Due  to  the  unknown 
normal  pressure  gradient,  the  reduction  of  the  measured  pitot 
pressures  requires  some  approximation. 

In  order  to  estimate  the  relative  importance  of  the 
terms  in  the  equations  of  motion  in  a  region  of  interest, 
one  introduces  the  following  nondimensional  variables. 

They  are  all  of  order  of  unity. 


X  =  x/X  , 

y  =  y/<s 

U  =  u/uro 

p  =  p/pto 

(19) 

P  =  P/PooU»  ' 

P  =  P/Pco 

Then,  from  the  continuity  equation  of  compressible 
flow,  one  knows  that  the  nondimensional  velocity  v  has 
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Figure  14.  Typical  Pitot  Surveys  of  Approaching  Boundary  Layer 


to  be  defined  as 


v  =  (v/uj  (X/6)  (20) 


With  the  nondiinensional  variables  defined  in 
Equations  (19)  and  (20),  the  nondiinensional  form  of 
the  Navier-Stokes  equations  for  the  steady  two-dimensional 
compressible  flow  are: 


__  9u _ 3u  3p  V'oo  X  2  3  __  3u 

P  u  —  +  P  v  —  =  -  —  +  -  (-)  —  (u  — ) 

3x  3y  3x  p^uJC  6  3y  3y 


^00  _  32u  P  32u  32v 

+  [U  -~yr  +  -  (—  +  ——) 

P^uJC  3x^  3  3  yT  3x3y 


2  3p  3u  3v  3y  3u  3y  3v 

- +  2  —  —  -t  ll  —  — )  (21) 

3  3x  3x  3y  3x  3x  3y  3x 


<S _ 3v _ 3v  X  3p 

-  [p  u  —  +  P  v  — ]  = - - 

X  3x  3y  5  3y 

X  32v  y  32u  32v 

+  — — —  {— )  [y  +  —  (■•„"_  +  “—y) 
P«,uooX  6  3y^  3  3x3y  3y 


2  3y  3u  3v  3v  Dy  3u  3y 

3  3y  3x  3y  3y  3y  3y  3x 


y^  6  32u  3v  3y 

+  -  (-)  [-35-  +  —  — ]  (22) 

P~.u™x  X  3x  3x  3x 

00  00 
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If  one  is  interested  in  the  development  of  the 
boundary  layer  along  a  whole  body,  then  X  is  the  body 
length,  6/L  <<  1;  and  Equations  (21)  and  (22)  can  be 
approximated  by: 


3u  3u  3p  y  L  2  3  _  3u 

P  u  —  +  P  v  —  - +  -  (-)  —  (y  — )  (23) 

3x  3y  3x  pjij,  5  3y  3y 


3p  6  2 

—  -  O(-)  (24) 

3y  L 


since  for  laminar  boundary  layer  flow 


6  2  1 
(-> - 

L  Re«,L 


Equations  (23)  and  (24)  reduce  to  the  boundary  layer 
equations . 

However,  if  one  is.  interested  in  a  region  whose 
axial  distance  is  comparable  to  the  boundary  layer  thick¬ 
ness,  then  X  5  6,  6/X  ~  1;  and  Equations  (21)  and  (22) 
can  be  approximated  by: 


3u  3u  3p  yw 

p  u  — -  +  pi  v  —  *  -  ■ —  +  -  (viscous  terms)  (25) 

3x  3y  3x  P.u.6 


3v  3v  3p  yw 

pu  —  +  p  v  —  -  -  —  +  - —  (viscous  terms)  (26) 

3x  3y  3y  P^u^fi 


since  for  laminar  flow 


1 


Pcc 


5 


then  Equations  (25)  and  (26)  become 


3u  _  3u  3p  6 

■p  u  +  p  v  —  =  -  —  +  —  (viscous  terms) 

3x  3y  3x  L 

3v _ 3v  3p  6 

p  u  —  +  pV“=-  —  +  —  (viscous  terms) 
3x  3y  3y  L 


(27) 


(28) 


Based  on  pitot  measurements  close  to  the  corner, 
one  finds  that  the  static  pressure  reaches  the  free  stream 
value  at  about  26  in  the  y  direction.  This,  and  the  fact 
that  the  boundary  layer  begins  to  accelerate  56  ahead  of 
the  corner,  suggests  that  Equations  (27)  and  (28)  are  the 
approximate  equations  of  motion  in  the  region  of  upstream 
influence. 

It  follows  from  Equations  (27)  and  (28)  that  the 
viscous  forces  on  a  fluid  element  are  of  the  order  6/L 
smaller  than  the  pressure  and  the  inertial  forces  in  the 
major  portion  of  the  expansion  (except  in  the  thin  viscous 
layer  adjacent  to  the  wall  where  non-slip  conditions  must 
be  satisfied).  Also,  the  momentum  equation  in  the  y 
direction  is  as  important  as  that  in  the  x  direction. 
Therefore,  it  will  be  a  good  approximation  to  treat  the 
acceleration  of  the  approaching  boundary  layer  and  the 
following  rapid  expansion  downstream  of  the  corner 
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isentropically.  Since  there  is  strong  entropy  gradient 
in  the  boundary  layer  before  acceleration,  the  stagnation 
pressure  is  constant  only  along  the  streamline. 

4.  APPROACHING  BOUNDARY  LAYER  ACCELERATION  ESTIMATED 
BY  STREAMTUBE  APPROXIMATION 

The  measured  pitot  pressure  profiles  are  analyzed 
based  on  the  isentropic  streamtube  approximation,  which 
assumes  that  total  pressure  is  conserved  along  each 
streamtube ,  i . e . , 

pox(13>  -  Pox(°3>  (29> 

i.u  i.u 

where  i  denotes  the  i  station  and  j  denotes  the  j 
streamtube.  The  notations  are  illustrated  in  Figure  15. 

i.U 

P  (o j )  is  the  stagnation  pressure  of  the  j  streamtube 

UA 

at  the  beginning  of  acceleration,  which  is  computed  from 

the  measured  pitot  pressure  under  the  boundary  layer  assumption. 

The  height  of  the  j*"*1  streamtube  from  the  wall  at 
i  station  is  computed  by 

j 

y (ij)  =  l  A (ij )  (30) 

j=0 

and  A(ij)  is  related  to  A(oj)  by  the  following  one¬ 
dimensional  isentropic  streamtube-area  relation  [41]. 


44 


45 


igure  15*  Isentropic  Streamtube  Approximation. 
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(31) 


where  A(oj)  is  the  initial  thickness  of  the  jth  streamtube. 

The  Mach  number  M(ij)  is  related  to  the  measured 

pitot  pressure  p  (ij)  by  the  normal  shock  equation  [41). 
oy 


Ppy^) 

Pox(i^ 


M2  (ij) 


Y 

Y-l 


Y-l  2 

1  +  -  M^(ij) 
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2Y 

Y+l 


M2  (ij) 


(32) 


Since  the  total  pressure  is  constant  along  each  streamtube, 
the  loci  of  constant  pitot  measurements  in  the  subsonic 
portion  of  the  boundary  layer  are  streamlines  (the  meas¬ 
ured  pitot  pressure  is  equal  to  the  stagnation  pressure) . 
The  physical  location  of  the  first  streamline  y(io)  is 
identified  somewhere  in  the  subsonic  portion  of  the  bound¬ 
ary  layer  from  the  pitot  pressure  measurements. 

Equations  (29)  through  (32)  can  be  combined  to  solve 
for  y(ij),  A(ij)  and  M(ij).  The  normal  pressure  gradient 
within  the  subsonic  portion  of  the  boundary  layer  was 
assumed  to  be  negligibly  small. 


The  Mach  number  profiles/  streamlines  and  normal 
static  pressure  distribution  at  stations  very  close  to  the 
separation  corner  were  calculated.  The  results  are  shown 
in  Figure  16.  Note  that  the  major  acceleration  takes 
place  at  the  lower  half  portion  of  the  approaching  bound¬ 
ary  layer.  The  pressure  gradient  normal  to  the  wall  be¬ 
comes  significant  in  the  supersonic  portion  of  the  bound¬ 
ary  layer  as  the  corner  is  approached.  The  static  pressure 
variation  in  the  supersonic  portion  of  the  boundary  layer 
is  found  to  be  approximately  linear. 

The  static  pressure  does  not  reach  the  free  stream 
value  at  the  outer  edge  of  the  approaching  boundary  layer, 
which  indicates  that  the  upstream  influence  was  not  con¬ 
fined  to  the  boundary  layer  but  appreciably  alters  the 
external  flow  as  well. 

It  is  interesting  to  observe  that  the  streamwise  and 
normal  pressure  gradients  near  the  corner  are  about  the 
same  magnitudes.  In  other  words,  (3p/3y)/ (3p/3x)  is 
approximately  unity.  It  does  not  indicate,  however, 
that  the  fluid  element  initially  close  to  the  surface 
will  accelerate  to  the  sonic  speed  at  the  corner,  nor 
does  it  show  any  large  curvature  of  the  streamlines  as 
predicted  from  the  surface  pressure  measurements.  This 
inconsistency  may  be  attributed  to  the  wall  effects  and 
the  local  probe  disturbance.  Also,  very  near  the  corner 
the  upstream  influence  to  the  probe  may  be  significant. 
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Static 


Humber  and  Static  Pressure  Distribution 


The  estimation  (as  shown  in  Figure  16)  not  only  gave  some 
indications  about  the  flow  field  near  the  corner,  but  it 
also  provides  the  initial  conditions  for  the  calculation 
by  method  of  rotational  characteristic  which  is  discussed 
in  the  following  chapter. 

The  difference  in  Mach  number  profile,  reduced  with 
and  without  static  pressure  gradient,  is  greater  at  the 
outer  edge  of  the  boundary  layer.  Also,  without  the 
static  pressure  correction,  the  reduced  total  pressure 
at  the  outer  edge  of  the  boundary  layer  is  greater  than 
the  free  stream  total  pressure.  Similar  estimation  for 
=  3.5,  Re^  L  =  1.0  x  10**  is  presented  in  Figure  17. 
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Static  Pressure  Distribution  at  M 


CHAPTER  V 


FLOW  DOWNSTREAM  OF  SEPARATION 

1.  PITOT  PRESSURE  MEASUREMENTS 

The  basic  features  of  the  entire  wake  of  a  back 
step  are  illustrated  in  Figure  18  by  the  pitot  pressure 
survey  for  Mm  =  3.5,  Re^  ^  =  1.0  x  10  ,  and  H  =  0.75 
inch.  The  lip  shock,  wake  shock  and  slip  line  can  be 
clearly  se.en,  and,  their  locations  agree  very  well  with 
the  shadow-graph  picture. 

The  pitot  profiles  in  the  lower  edge  of  the  free 
shear  layer  and  recirculating  region  are  shown  in  Figure 
19.  Note,  the  existence  of  the  low  pressure  which  dies 
away  downstream  along  a  line  which  is  suspected  to  be  the 
dividing  streamline  from  hot-wire  resistance  measurements.* 
This  may  be  the  result  of  viscous  effect  on  the  pitot  probe 
at  very  low  Reynolds  number;  and  it  is  not  likely  that 
this  low  pressure  is  due  to  the  probe  interference  or 
the  difference  in  flow  direction  with  respect  to  the  probe. 

The  pitot  pressure  distributions  in  the  immediate 
neighborhood  of  the  separation  corner  are  presented  in 
Figure  20,  where  the  lip  shock  and  the  initial  develop¬ 
ment  of  the  free  shear  layer  are  evident.  Since  the 

*This  point  is  discussed  further  in  Chapter  IX. 
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Backstep 


stagnation  pressure  is  relatively  lower  in  the  shear  layer 
expanded  from  the  lower  portion  of  the  boundary  layer  at 
the  corner,  the  measured  pitot  pressure  variation  acror-s 
the  lip  shock  will  be  very  small  if  the  lip  shock  origi¬ 
nates  at  the  corner  and  passes  through  this  region. 

Therefore,  it  is  rather  difficult  to  detect  the  beginning 
of  tAe  lip  shock  by  pitot  pressure  surveys.  A  very  dense 
pitot  pressure  measurement  about  one  boundary  layer  thick¬ 
ness  downstream  of  the  corner  appears  to  support  this 
reasoning  (as  shown  in  Figure  21) . 

Since  there  are  entropy  gradients  and  streamwise  and 
lateral  pressure  gradients  in  the  expanding  boundary  layer, 
the  pitot  measurement  in  this  region  does  not  provide 
sufficient  information  for  quantitative  analysis.  To 
supplement  the  pitot  measurement,  hot-wire  measurements 
were  attempted  in  this  region.  Unfortunately,  in  the  regions 
of  higher  velocity  near  the  corner  and  at  reattachment, 
severe  wire  breakage  was  encounted.  For  this  reason,  no 
hot  wire  results  were  obtained  in  these  regions,  and  only 
the  cavity  flow  region  was  explored  in  this  way.  These 
results  are  discussed  in  Chapter  VIII.  However,  comparisons 
of  the  measured  pitot  profiles  with  the  results  of  theoretical 
calculations  may  provide  some  valuable  information. 

2.  COMPARISONS  WITH  RESULTS  OF  PRANDTL-MEYER 

EXPANSION 

The  pitot  pressure  measurements  have  been  compared 
with  calculations  from  the  Prandtl-Meyer  equation  for  the 
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following  cases: 


M 

Re  . 

00 

00, L 

2.5 

1.06  x  10 

3.5 

1.06  x  10 

5.07 

1.9  x  106 

Since  the  Prandtl-Meyer  theory  is  restricted  to  irrotational 
supersonic  flow,  it  is  expected  that  the  existence  of  the 
boundary-layer  will  cause  departures  from  the  simple  theory. 
The  comparisons  are  shown  in  Figs.  22,  23,  24, 
respectively. 

Scherberg  and  Smith  [18,  25]  and  Hama  [8]  observed 
the  phenomenon  of  an  upstream  influence  of  the  lower  base 
pressure  on  the  approaching  boundary-layer  before  the  expan¬ 
sion  corner.  This  point  has  been  confirmed  by  this  study, 
which  reports  in  greater  detail  the  pressure  measurements 
in  this  region.  A  correlation  of  the  surface  pressure  dis¬ 
tribution  upstream  of  the  expansion  corner  has  been  dis¬ 
cussed  in  the  previous  chapter.  In  addition,  the  over¬ 
expansion  of  the  external  inviscid  stream  ahead  of  the  lip 
shock  has  been  clearly  demonstrated. 

The  lip  shock  is  embedded  in  the  expanded  boundary- 
layer  close  to  the  corner  and  penetrates  into  the  region 
of  vanishing  vorticity.  The  pitot  pressures  at  the  end  of 
the  expansion  do  not  achieve  the  values  predicted  by  the 
Prandtl-Meyer  theory,  but  rather  approach  it  after  about 
10  to  15  boundary-layer  thickness  downstream.  Also  beyond 
this  station  (see  Figs.  22-24)  there  is  a  region  of  nearly 
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Figure  22.  Comparison  of  Measured  Pitot  Profiles  with  Prandtl-Meyer 
Expansion  at  M  =2.5. 


X  =  O  l6  inch 


Comparison  with  Prandtl-Meyer  Expansion  at 


uniform  inviscid  flow  immediately  downstream  of  the  lip 

shock  as  indicated  by  the  constant  pitot  pressures,  which 

are  nearly  equal  to  the  Prandtl-Meyer  value.  Scherberg 

and  Smith  [18]  have  pointed  out  in  their  paper  that  a 

constant  pitot  pressure  reading  does  not  necessarily 

indicate  a  uniform  flow.  Thus,  they  indicate  that  a 

sensibly  constant  total  head  reading  ahead  of  an  oblique 

shock  can  admit  changes  in  the  component  velocity  and 

static  pressures  such  that  the  shock  strength  need  not 

be  constant.  Therefore,  a  careful  check  was  made  by 

examining  a  series  of  shadow  graphs  with  Kama's  method. 

Within  the  accuracy  of  Kama's  method,  the  Mach  number 

behind  the  lip  shock  is  found  to  be  very  nearly  a  constant 

value  at  different  stations  for  a  given  free  stream  Mach 

number  and  Reynolds  number  (but  the  flow  direction  need  not 

be  constant  [18]).  This  is  more  noticable  for  the  lower 

free  stream  Mach  number  cases  than  for  the  higher  free 

stream  Mach  number  case  (such  as  M  =5.0).  Since  the  static 

00 

pressure  in  this  region  is  very  close  to  the  base  pressure, 
it  is  clear  that  the  external  inviscid  flow,  upstream  of 
the  lip  shock,  overexpands  to  a  pressure  lower  than  the 
base  pressure.  Except  close  to  the  corner,  the  lip  shock 
is  at  a  considerable  distance  above  the  outer  edge  of  the 
free  shear  layer. 

It  is  obvious  and  expected  that  the  Prandtl-Meyer 
theory  fails  to  describe  the  flow  field  properly  in  such 
cases  because  of  the  boundary-layer  shed  from  the  body 
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surface.  To  accurately  describe  the  flow  field,  the  in¬ 
coming  boundary-layer  cannot  be  neglected  in  the  expansion- 
separation  process  (more  examples  will  be  given  later) . 
Furthermore,  a  similar  phenomenon,  i.e.,  a  dip  in  the  sur¬ 
face  pressure  before  a  sharp  corner  in  a  high  subsonic 
flow  was  reported  by  Rubin  [43]. 

The  expansion  of  the  supersonic  boundary-layer 
around  a  separation  corner  has  been  treated  with  inviscid 
rotational  theory;  and  the  cause  of  lip  shock  formation 
was  described  as  the  relatively  slow  coalescence  of  the 
reflected  wavelets  [42].  It  is  of  interest  to  compare  the 
pitot  measurements  with  the  results  of  an  inviscid  rotational 
method  of  characteristics  calculation;  and  this  motivated 
the  following  investigation. 

3.  COMPARISONS  WITH  RESULTS  OF  ROTATIONAL  CHARACTERISTICS 
CALCULATION 

The  pitot  pressures  computed  by  the  method  of  rotational 
characteristics  were  provided  by  Matz  [44]  and  are  shown 
in  Figs.  25  and  26  for  Mro  =  3.5  and  Mm  =  5.07  respectively. 

The  initial  conditions  for  the  computation  at  the  separation 
corner  were  estimated  [26,  45]  based  on  the  measurement  of 
the  approaching  boundary-layer  before  the  corner.  The  super¬ 
sonic  portion  of  the  boundary-layer  at  the  corner  was  expanded 
to  the  base  pressure  (either  a  measured  value  or  an  estimated 
one) ,  while  the  expansion  of  the  inner  boundary-layer  at 
subsonic  speed  is  neglected.  The  base  pressure  has  been 
studied  in  detail,  as  outlined  in  Chapter  VI,  and 


permits  us  to  estimate  the  base  pressure  fairly  accurately. 

Very  good  agreement  was  found  in  the  expanding 
* 

boundary  layer  for  both  cases.  As  shown  in  Fig.  25,  the 
pitot  pressure  profiles,  the  computed  lip  shock  positions 
and  the  pitot  pressure  changes  across  the  lip  shock  agree 
remarkably  well  with  the  measurements.  This  appears  to 
support  the  contention  that  the  formation  of  the  lip  shock 
is  an  inviscid  rotational  field  phenomena-as  originally  pro¬ 
posed  by  Charwat  and  Yakura  [24]  and  investigated  by  Weinbaum 
[42],  The  effect  of  viscosity  appears  to  be  negligibly 
small  compared  to  the  change  in  pressure  during  the  rapid 
expansion  of  the  supersonic  boundary-layer  downstream 
of  a  separation  corner. 

Typical  profiles  of  Mach  number,  static,  total  and 

pitot  pressures  obtained  from  the  rotational  characteristic 

calculations  for  two  stations  are  shown  in  Fig.  27  for  M  = 

00 

£ 

3.5,  Re  .=1.0x10  and  H  =  0.75  inch.  The  reduced  Mach 

00, 

number  from  hot-wire  measurements  are  also  included  for 
comparison.  Note  the  lower  static  pressures  and  the  corres¬ 
ponding  overexpansion  (compared  to  the  Prandtl-Meyer  expansion 
to  the  same  base  pressure)  just  before  the  lip  shock.  At 
x  =  0.8  inch  (roughly  a  step-height  distance  downstream  of 
the  corner) ,  the  higher  total  pressure  at  the  lip  shock 
makes  the  jump  in  pitot  pressure  across  the  lip  shock  larger 
and  much  more  observable  than  that  at  x  =  0.4  inch  station; 
even  though  the  Mach  numbers  and  static  pressures  before  and 

after  the  lip  shock  are  about  the  same  for  both  stations. 

*It  is  defined  herein  that  the  "shear-layer"  is  the  "expanded 
boundary  layer"  downstream  of  the  corner. 
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Comparison  of  Measured  and  Computed  Pitot  Pressures 


Typical  Results 


The  calculated  pitot  pressure  normalized  with  the 
total  pressure  V  /V.,  and  the  Mach  number  immediately 

OX  v 

upstream  (denoted  by  sub  x)  of  the  lip  shock,  M  ,  are 

X 

shown  in  Fig.  28  and  indicate  that  within  a  distance  of 

one  step  height,  the  lip  shock  is  embedded  in  the  vorticity 

layer  originated  from  the  expanded  boundary-layer,  and 

that  the  Mach  number  variation  along  the  lip  shock  is  small. 

The  static  pressure  jump  P  /P  ,  across  the  lip  shock  is 

sy  sx 

also  shown  (subscript  y  indicates  conditions  downstream 
of  the  shock) . 

4.  COMPARISON  OF  MEASURED  PITOT  PROFILES  AT  DIFFERENT 
FLOW  CONDITIONS 

The  pitot  pressure  measurements  at  =  3.5  and 

Re  .  =  1.0  x  106  and  1.72  x  10®  are  shown  in  Fig.  29, 

00,1, 

where  the  boundary-layers  are  laminar  before  separation  in 
both  cases.  The  pitot  pressure  is  again  normalized  by  its 
corresponding  free  stream  total  pressure.  The  flows  ex¬ 
panding  from  the  external  inviscid  streams  under  these  two 
conditions  are  very  nearly  the  same  (see  x  =  0.1  inch 
station  in  Fig.  29)  and  they  approach  asymptotically  to 
the  flow  field  predicted  by  the  Prandtl-Meyer  theory. 
However,  the  flow  turns  through  an  angle  that  is  consistent 
with  its  corresponding  base  pressure,  the  direction  and 
profile  of  the  free  shear  layers  are  very  different  down¬ 
stream.  As  a  reference  and  a  comparison,  the  value  cal¬ 
culated  by  Prandtl-Meyer  expansion  to  the  base  pressure 

g 

corresponding  to  that  of  Rew  L  =  1.0  x  10  'is  included  in 
Fig.  29. 
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Comparisons  of  Pitot  Pressure  Profiles  for  Differen 
Reynolds  Numbers  and  Inviscid  Flow  Theory  at  M  =  3 


Fig.  30  shows  the  measurements  for  the  M  =2.5, 

CO 

Re  .  =  1.0  x  10  and  1.82  x  10°  flows.  The  considerable 

00, l 

deviation  in  the  expanding  inviscid  streams  are  probably 
due  to  the  fact  that  for  Re^  ^  =  1.82  x  10°,  the  boundary- 
layer  is  transitional  near  the  corner,  and  the  boundary 
layer  thickness  increases  as  the  Reynolds  number  increases. 

For  a  given  free  stream  Mach  number,  the  lip  shock  is 
directed  downward  at  the  high  Reynolds  number  but  turns 
upward  and  decreases  in  strength  as  the  Reynolds  number 
decreases.  The  reduced  Mach  number  profiles,  corresponding 
to  the  flow  conditions  of  Fig.  30,  is  shown  in  Fig.  31.  The 
calculation  was  based  upon  the  assumptions  that:  (i)  above 
the  expanded  boundary-layer  the  total  pressure  is  equal  to 
the  free  stream  total  pressure,  (ii)  underneath  the  lip 
shock  the  static  pressure  .is  equal  to  the  base  pressure. 

It  is  quite  evident  from  these  observations  that  the 
post-expansion  flow  field  is  strongly  dependent  upon  the 
nature  of  approaching  boundary-layer  at  the  corner  and  upon 
the  base  pressure,  and,  both  are  a  function  of  Re^  L  and 
other  factors. 

5.  ON  FREE  STREAM  MACH  NUMBER  EFFECT 

The  measured  pitot  pressure  profiles,  one  inch  down¬ 
stream  of  the  step,  are  shown  in  Fig.  32  for  Re^  ^  =  1.0  x 
6 

10  and  free  stream  Mach  number  of  2.5,  3.5  and  5.0  respectively. 
The  measurements  are  normalized  by  the  free  stream  stagnation 
pressure.  These  profiles  are  presented  in  different  scales. 

The  upstream  influence  in  the  outer  inviscid  expansion  and 
the  departure  of  the  separated  profiles  from  the  Prandtl-Meyer 
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prediction  are  more  severe  at  the  high  free  stream  Mach 
number  flow  as  can  be  seen.  The  lip  shock  also  turns 
outward  as  the  free  stream  Mach  number  increases.  Note 
also  that  the  thickness  of  the  free  shear-layer  increase 
as  the  Mach  number  increases  while  the  region  of  the  uni¬ 
form  inviscid  flow  downstream  of  the  lip  shock  decreases 
(Fig.  33  and  also  see  Fig.  22,  23,  24).  At  =  5.0, 
within  a  distance  of  one  step  height  downstream,  there  is 
no  evidence  of  the  existence  of  a  uniform-flow  region  be¬ 
tween  the  shear-layer  and  the  lip  shock.  The  lip  shock  is 
embedded  in  the  shear-layer  for  a  greater  distance. 

6.  INFLUENCE  OF  STEP  HEIGHTS 

The  effects  of  step  height  on  the  expanded  and 
separated  flow  field  are  shown  in  Fig.  33  for  M^  =  3.5, 

Re  .  =  1.05  x  10°.  The  step  heights  are  0.75  inch  and 

00,  lj 

0.44  inch  respectively.  It  seems  that  the  upstream  influence 
due  to  different  step  heights  is  negligible  in  that  the 
expanded  profiles  above  the  lip  shock  are  very  similar 
initially.  The  location  of  the  lip  shocks  are  strongly 
related  to  the  base  pressures.  The  differences  in  the 
separated  boundary-layer  are  mainly  due  to  the  degree  of 
expansion.  Also,  the  results  of  rotational  characteristics 
calculation  are  included,  the  initial  conditions  at  the 
separation  corner  are  assumed  identical  to  that  of  Fig.  26. 
This  provides  further  evidence  that  the  rotational  character¬ 
istics  are  adequate  to  predict  the  expanded  flow  field. 
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CHAPTER  VI 


A  NEW  BASE  PRESSURE  CORRELATION 

As  outlined  in  Chapter  II,  detailed  pressure  distri¬ 
butions  on  the  step  were  measured.  This  base  pressure 
(measured  on  the  reattachment  wall  close  to  the  base  of 
the  step)  is  a  function  of  Mach  number  and  Reynolds 
number.  A  correlation  of  this  data  has  been  attempted  as 
outlined  below. 

In  a  local  separation  region,  a  Reynolds  number  based 
on  the  characteristic  length  of  the  body  is  not  necessarily 
related  to  the  physical  aspects  of  the  local  conditions-. 

For  a  base  flow,  it  is  easily  sensed  that  the  characteristic 
lengths  are  the  shear  length  S-m  and  the  mixing  width  Am  in 
the  recompression  region  (see  Fig.  34).  An  attempt  has 
been  made  to  correlate  the  base  pressure  with  these  two 
lengths,  and  a  very  encouraging  correlation  was  obtained. 

The  only  difficulty  is  that  the  lengths  -tm  and  the  Am  are 
not  known  a  priori.  Therefore,  these  two  lengths  must  be 
related  to  the  known  body  lengths;  namely,  the  model  length 
upstream  of  a  two  dimensional  rearward  facing  step  L  and 
the  step  height  H.  Physically,  the  length  L  is  related  to 
the  boundary  layer  thickness  5  prior  to  separation,  and 
this  parameter  will  influence  the  mixing  width  A  .  The 
step  height  H  is  related  to  the  shear  layer  length  tm, 
when  the  other  flow  parameters  are  held  constant. 
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om  this 


The  ratio  of  inertia  to  viscous  forces  in  the  shear 
layer,  or  the  physically  correct  reduced  Reynolds  number 
Rec  in  this  region,  may  be  expressed  as: 
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where  and  are  taken  as  the  characteristic  dimensions 
of  the  mixing  region  as  indicated  above.  The  estimation  of 
^m  and  \  as  °^ten  difficult.  If  one  assumes  that 
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where  the  exponent  n  is  a  function  of  free  stream  Reynolds 

number  based  on  L  and  the  free  stream  Mach  number  M  .  The 

00 

exponent  n  may  depend  strongly  on  whether  the  shear  layer 

is  laminar,  transitional  or  turbulent.  Based  on  equations 

(34)  and  (35),  the  reduced  Reynolds  number  Re  of  Equation 

c 

(33)  becomes, 
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If  one  also  neglects  the  dependence  of  6  on  Mach 


number,  then, 
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and 


u  m  +  2n  -  2  T  m  -  1  (38) 

«*o  -  <  I  >  (h> 

where  for  very  small  pressure  gradient  flow,  m  is  (approxi¬ 
mately)  2  or  5  depending  on  whether  the  boundar  layer  on 
the  body  is  laminar  or  turbulent  respectively.  The  reduced 
Reynolds  number  is  thus  related  to  6,  H  and  L,  which  are 
the  important  parameters  of  the  base  flow.  This  reduced 
Reynolds  number  is  very  likely  to  be  a  correlation  para¬ 
meter  for  the  base  pressure. 

Based  on  their  mixing  theory,  Crocco  and  Lees  [4]  were 
able  to  show  the  important  trends  in  the  variation  of  base 
pressure  with  Reynolds  number,  L/H  and  the  nature  of  the 
boundary  layer.  They  argued  that  the  base  pressure  is  a 
function  of  the  mixing  rate  and  of  5.  The  mixing  rate  and  6 
depend  not  only  on  Reynolds  number  but  also  on  the  location 
of  transition  (whether  this  be  on  the  body  or  in  the  separated 
region).  For  an  intermediate  wake  transition  region,  the 
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rate  of  upstream  movement  of  transition  in  the  wake  with 
increasing  Reynolds  number  is  governed  largely  by  the 
parameter  L/H.  For  fully  laminar  or  turbulent  flow,  where 
L/H  is  not  as  important,  the  base  pressure  ratio  then 
depends  mainly  on  the  ratio  6/H;  as  correlated  by  Chapman 
[45],  For  transitional  wake  flow,  m  =  2  is  appropriate, 
and  if  it  is  assumed  that  n  =  0.&*,  the  correlation  parameter 
reads 
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Because  of  the  contrast  in  dependence  of  the  base  pressure 
on  the  mixing  rate  and  on  the  boundary  layer  thickness  (as 
pointed  out  by  Crocco  and  Lees  [4]),  it  is  proposed  to  nor¬ 
malize  the  base  pressure  ratio  with  the  correlation  para¬ 
meter.  The  correlation  is  shown  in  Fig.  35  for  the  data 
the  range  of  the  data  is  as  follows  (see  Table  1): 

0.05  x  106  <  Re  .  <  1.8  x  106 

00,  Jj 

2.0  <  M  <3.55 

CO 

5.33  <  L/H  <  31.4 

0.0072  inch  <  H  <  0.75  inch 
The  proposed  base  pressure  correlation  for  a  two 
dimensional  rearward  facing  step  is, 

♦This  value  has  been  found  empirically. 
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Table  1.  Base  Pressure  Correlation. 


Symbol 

M 

00 

H,  inch 

L/H 

Investigators 

© 

2.50 

0.7500 

5.33 

Present 

A 

3.50 

0.7500 

5.33 

Present 

o 

2.50 

0.4400 

9.03 

Smith 

A 

3.50 

0.4400 

9.03 

Smith 

<£> 

2.00 

0.0100 

20.00 

Chapman 

O 

2.00 

0.1000 

20.00 

Chapman 

X 

2.25 

0.0072 

31.40 

Rom 

SJ 

3.55 

0.0072 

31.40 

Rom 

0 

2.25 

0.0165 

13.20 

Rom 

0 

3.55 

0.0165 

13.20 

Rom 

0 

2.25 

0.0256 

13.00 

Rom 

0 

3.55 

0.0256 

13.00 

Rom 

0 

2.25 

0.0335 

7.40 

Rom 

0 

3.55 

0.0335 

7.40 

Rom 

0 

2.25 

0.0790 

8.00 

Rom 

0 

3.50 

0.0790 

8.00 

Rom 

0 

2.25 

0.0790 

15.00 

Rom 

0 

3.50 

0.0790 

15.00 

Rom 

0 

5.00 

0.7500 

5  .00 

Present 

Source:  H.  E.  Smith,  "The  Plow  Field  and  Heat 
Transfer  Downstream  of  a  Rearward  Pacing  Step  in  Super' 
sonic  Flow,"  Aerospace  Research  Laboratory  Report 
No.  67-0056  (Wright-Patterson  Air  Force  Base,  Ohio, 
March,  1967) . 

Source:  D.  R.  Chapman,  D.  M.  Kuehn,  and 
H.  K.  Larson,  "Investigation  of  Separated  Flows  in 
Supersonic  and  Subsonic  Streams  with  Emphasis  on  the 
Effect  of  Transition,"  National  Advisory  Committee  for 
Aeronautics  Report  No.  1356  (Washington,  D.  C.,  June, 
1958) . 


Source:  J.  Rom,  "Supersonic  Flow  over  Two- 
Dimensional  and  Axially  Symmetric  Backward  Facing.  Steps," 
TAE  Report  No.  33  (Department  of  Aeronautical  Engineering, 
Technion-Israel  Institute  of  Technology,  Haifa,  Israel, 
March,  1964). 
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(40) 


by  o 

p-  =  Re  exp(-  6.97  +  0.745D  -.0.1136D  ) 

00 

where  D  =  -tn(Rec) 

The  comparison  of  measured  base  pressure  ratio  and 
the  base  pressure  ratio  computed  by  Equation  (40)  is  shown 
in  Fig.  36,  which  indicates  that  the  correlation  is  indeed 
a  reasonable  one. 

Similar  correlation  for  a  6  degree  wedge  [8]  is 
shown  in  Fig.  37.  Note  that  in  both  Figs.  35  and  37,  the 
body  boundary  layer  is  laminar  before  separation,  and  that 
this  correlation  parameter  applies  to  fully  laminar  flow 
as  well. 
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Pjj/Po,  (Computed  from  Correlation) 


Figure  36.  Comparison  of  Measured  Base  Pressure  and 
Base  Pressure  Computed  from  Correlation. 


CHAPTER  VII 


THE  INFLUENCE  OF  AN  OVERHANG 
MOUNTED  AT  THE  SEPARATION  CORNER 

It  has  been  postulated  that  the  dividing  streamline 
actually  turns  sharply  around  the  separation  corner  and 
then  separates  from  the  surface  of  the  step  [8,  57^  58], 
Whatever  the  mechanism  of  the  separation,  it  is  possible 
to  change  the  conditions  at  separation  by  changing  the 
step  geometry  in  this  region.  In  order  to  determine  the 
effect  of  different  conditions  at  separation,  a  thin 
(0.02  inch)  overhang  of  various  lengths  was  added  to  the 
basic  model.  The  inset  on  Fig.  39  shows  the  overhang 
attachment. 

A  theoretical  model  that  assumes  a  ’dead-air*  region 
downstream  of  the  step  would  not  differentiate  between 
the  plain  step  and  the  step  with  overhang.  Indeed,  for 
low  supersonic  flow  Nash  et  al  [55]  reported  no  measurable 
effect  on  base  pressure  by  the  introduction  of  an  overhang 
configuration.  The  same  comment  was  not  true  in  subsonic 
flow;  but  a  discussion  of  that  flow  is  not  the  subject  of 
the  present  work.  If  there  are  any  differences  between 
the  two  flows  it  must  result  from  changes  in  the  geometry 
of  the  streamlines  in  the  "dead-air"  region  near  the  se¬ 
paration  point.  This  could  have  some  effect  on  the  mixing 
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in  the  shear  layer  due  tc  modifications  of  the  velocity  profiles. 

It  is  considered  that  this  effect  will  not  be  too  large  [59]. 

The  results  shown  in  Figures  38  and  39  are  not  immedi¬ 
ately  interpretable.  Thus  at  a  Mach  number  of  3.  5  (Figure  38), 
there  is  no  noticeable  change  in  base  pressure  with  configura¬ 
tion.  The  base  pressure  correlation  (eauation  40)  shows  a  very 
small  change  in  base  pressure  over  the  Reynolds  number  of  the 
test  data.  For  a  Mach  number  of  2.  5  (Figure  39),  the  data  shows 
an  increase  of  base  pressure  with  the  introduction  of  an  overhang; 
which  is  opposite  to  the  effect  anticipated  from  the  change  of 
Reynolds  number.  This  inconsistency  of  the  data  with  change  in 
Mach  number  suggests  that  further  studies  should  be  made. 
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CHAPTER  VIII 


HOT-WIRE  MEASUREMENTS  IN  RECIRCULATING 
REGION  AND  FREE  SHEAR  LAYER 

1.  GENERAL  DISCUSSION 

Hot-wire  anemometer  probes  were  used  to  explore  the 
basic  structure  of  the  recirculating  region?  The  major 
difficulties  encountered  in  these  measurements  were: 

1.  The  Reynolds  number  based  on  the  wire  diameter 
was  very  low  (about  0.2  to  9). 

2.  Abrupt  velocity  and  total  temperature  gradients 
were  present  in  the  free  shear  layer. 

Considerable  effort  has  been  expended  in  developing  a 
satisfactory  data  reduction  technique/  which  was  presented 
in  Chapter  III. 

The  continuous  hot-wire  resistance  traces  at  con- 
-3 

stant  current  (7  *  10  amp)  are  presented  in  Figure  41. 
The  striking  changes  in  resistance  are  due  to  the  combined 
effects  of  changes  in  heat  transfer  rate  and  the  adiabatic 
wire  temperature  as  indicated  by  Equation  (12).  However, 
calculation  based  on  Equation  (13)  shows  that  at  small 
current  the  wire  measures  essentially  the  adiabatic  wire 
temperature,  which  is  the  product  of  recovery  factor  n  and 
the  local  total  temperature  TQ,  i.e.,  Taw  =  nTQ. 


^Analytical  work  was  performed  by  Scherberg  (15)  who  found  the  core 
streamline  structure  with  essentially  zero  pressure  gradients  in  the 
y  direction. 
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Dividing  Streamline 


In  the  region  of  the  measurements,  the  recovery  factor 
of  the  wire  varies  from  1.0  to  about  1.08,  depending  on 
the  local  Mach  number  and  Knud sen  number.  The  local  total 
temperature  in  the  free  shear  layer  and  recirculating 
region  is  not  uniform.  Therefore,  the  measured  relative 
minimum  resistances  (minimum  Taw)  are  primarily  due  to 
either  low  total  temperature  or  lower  recovery  factor. 

The  minimum  resistance  further  away  from  the  wall,  as 
shown  in  Figure  40,  is  due  to  the  lower  TQ  along  the 
dividing  streamline.  This  point  will  be  discussed  in 
greater  detail  later.  While  the  other  minimum 
resistance  is  due  to  the  lower  recovery  factor  at  the 
lower  Mach  number  (or  minimum  velocity) ,  the  recovery 
factor  decreases  to  unity  regardless  of  Knudsen  number 
in  the  range  investigated  as  the  Mach  number  approaches 
zero  (Figure  59,  Appendix  F) . 

Figure  41  reports  the  typical  static  pressure, 
pitot  pressure  and  hot-wire  resistance  measured  at  x/H  = 
1.14.  The  comparison  shows  the  necessity  of  employing 
hot-wire  techniques  for  measurements  in  the  lower  portion  of 
the  free  shear  layer  and  in  the  recirculating  region. 

There  is  some  doubt  as  to  the  effect  of  turbulence 
on  the  wire  resistance  measurements,  since  part  of  the  free 
shear  layer  is  transitional.  However,  estimations  indi¬ 
cate  that  this  is  negligible,  and  this  is  discussed 
further  in  Appendix  F. 
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inch,  L  =  4.0  inches,  H  =  0.75  inch 


Figure  41.  Comparison  among  Measured  Hot-Wire,  Pitot  and  Static  Pressure  Profiles 


2.  REDUCED  RESULTS 


The  variation  of  various  flow  properties  in  the 
free  shear  layer  and  the  recirculating  region  are  shown 
in  Figure  42.  h  and  T  are  directly  reduced  from  hot¬ 
wire  measurements.  The  broken  line  represents  values 
reduced  from  continuous  resistance  traces  which  are  not  as 
accurate  as  the  discrete  measurements.  By  combining  h, 

T  „  and  p  ,  the  other  flow  parameters  such  as  Mach  number 
and  total  temperature  are  then  determined  by  the  method 
described  in  Chapter  III.  The  Nusselt  number  decreases 
very  sharply  in  the  free  shear  layer  and  it  varies  also 
in  the  recirculating  region.  The  Knudsen  number  increases 
from  0.2  at  the  outer  edge  of  the  free  shear  layer  to  1.8 
in  the  recirculating  region  indicating  that  the  hot-wire  was 
in  the  transitional  (slip  flow)  regime.  The  distribution  of 
Reynolds  number  based  on  wire  diameter  (Rewd  0(1) ,  which  is 
essentially  the  mass  flux  distribution,  is  also  shown  in  Fig.  42. 

The  wire  recovery  temperature  and  convective  heat 
transfer  rate  for  various  stations  in  the  recirculating 
region  and  shear  layer  are  shown  in  Figures  43  and  44 
respectively.  The  basic  features  of  the  recirculating 
region  can  be  seen  more  clearly  from  the  following  reduced 
results. 

The  Mach  number  mapping  is  presented  in  Fig.  45 
with  the  dividing  steamline  and  an  approximate  indication 
of  the  zero  "u"  component  of  velocity  shown.  Although  the 
Mach  number  is  positive  below  the  latter  line,  the 
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by  Point  Measurements 
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Figure  42.  Typical  Reduced  Results  from  Hot-Wire  Measurements 


culating  Region 


Number  in  Recirculating  Region  and  Shear  Layer 


flow  direction  is  reversed.  Note,  the  acceleration  of 
the  flow  along  the  dividing  streamline.  Since  the  recir¬ 
culating  flow  is  of  finite  velocity  (but  the  Mach  number  can 
be  as  high  as  0.8),  the  interaction  between  the  separated 
boundary  layer  and  the  recirculating  flow  could  be  impor¬ 
tant. 

The  total  temperature  mapping  is  shown  in  Figure 
46.  The  total  temperature  increases  rather  rapidly  from 
440°R  close  to  the  wall  to  500°R  within  one  step  height 
distance  downstream,  while  the  rest  of  the  recirculating 
region  is  at  a  nearly  constant  total  temperature  of  500°R 
(free  stream  total  temperature  is  560 °R) . 

The  development  of  the  temperature  profile  will  be 
discussed  further  in  Chapter  IX  from  a  theoretical  stand¬ 
point,  and  further  deliberations  will  be  delayed  until 
that  time. 
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Figure  46.  Stagnation  Temperature  in  Recirculating  Region 


CHAPTER  IX 


AN  EXTIMATION  OF  THE  SHEAR-LAYER 
TEMPERATURE  DISTRIBUTION 


If  the  shear  layer  is  ’thin'  compared  to  its  length, 
it  can  be  assumed  that  the  boundary-layer  equations  hold 
in  a  coordinate  system  aligned  with  the  shear-layer  (see 
Fig.  47),  then: 


S  (pu>  +  (pv)  -  0 


(41) 


=  a  /.,au 


P  fit  + 


Dt  dx 


Sr  (^> 


(42) 


D(TC  ) 
P  "Dt 


3  ( k  St 


(*£)♦«* 


3?'  S7 


+  Ud£ 
dx 


(43) 


Introducing  the  total  temperature: 


Tq  =  T  +  u  /2  Cp 


(44) 


The  energy  equation  then  becomes : 


r  ®  (T  r  >  -  —  ( kBTo )  ,  u  S  1  \  -Sul 

P  m  (ToV  "  3y  \-^rl  +  -ft  [  H*1  -  p-> 


1  V  ,du 


2  N  r,  2  SC  1 

+  ud  -  —)®  +  ^[2^72-  ^ 


+  57  2C 


(43a) 


101 


102 


Figure  47.  Coordinate  System  for  Shear  Layer  Calculations. 


where : 


is  the  Prandlt  Number. 


pr  =  ^ 

Restricting  attention  to  flows  for  which  p?  =  1  and 
assuming  variations  of  the  specific  heat  at  constant 
pressure  are  small,  it  follows  that: 


(45) 


Introducing  a  Kirchhoff  transformation  with  new 
dependent  variable  tq  (Tq) ,  i.e.: 


(46) 


Then,  equation  45  becomes: 


for  unit  Prandtl  Number. 

Note  that  for  zero  pressure  gradient;  the  momentum 
equation  reduces  to: 


„  Du  _  d  (\x  bu) 

p  Dt  =  57  I  57/ 


(42a) 


so  that 


to  =  A  +  B  u  (48) 

where  A  and  A  are  constant.  This  result  is  only  valid  if 
Cp  =  const.,  Pr  =  1  and  -2.  =  o.  In  general,  this  result 
can  only  be  applied  to  the  asymptotic  shear  layer;  since 
the  initial  velocity  and  temperature  profiles  will  not 
be  related  in  this  way.  In  the  developing  shear  layer  it 
is  necessary  to  specify  both  initial  velocity  and  tempera¬ 
ture  profiles  and  solve  the  equations  45  and  42a  separately. 
Since  both  equations  are  of  the  same  form  (assuming  Pr  =  1 
in  the  energy  equation),  the  solution  for  only  one  need 
be  discussed. 

1.  APPROXIMATE  SOLUTION  OF  THE  EQUATION  (42a) 

For  steady  flow  the  equation  reads: 

„  du  ,  ...  du  _  c>  /  ix  du  ) 

Pu  5*  +  Pv  dy  "  37  (  3y / 

Introducing  the  new  variables  (after  Stewartson  [56]): 

dy  (49) 

P  00 

X  =  x 


104 


Then  if  the  stream  function  ip  is  given  by: 


as  suggested  by  the  continuity  equation,  then  it  follows 
(with  the  additional  assumption  that  the  viscosity  co¬ 
efficient  is  proportional  to  the  temperature)  that: 


-yhJ  M  a?  a23” 

Sx  ~  cix  <^Y 


where: 

X  =  M 

*  By 


(50) 


Neglecting  the  second  term  in  equation  (50)  (which 
is  equivalent  to  the  assumption  (pv  ^)/(pu  )«  1)  and 

replacing  the  coefficient  '7*  in  the  first  term  by  a 
constant  mean  value  7,  there  results  a  diffusion  equation: 


*  ay  a2y 

5  ax  = 


(51) 


in  the  variables  X,  Y, 
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Equation  (51)  solves  to  give: 


^(X,  Y)  =  -± 


”5  (o,  t))  exp  (-02)d£ 


where : 


n  =  1  -  ^fv  ~ 

7 


In  this  study,  7  has  been  taken  to  be  the  value  in 
the  free  s'  ?eam,  and  the  development  of  the  profile  is 
given  in  terms  of  the  initial  profile  (o,y)  and  the 
distance  downstream.  In  the  case  of  a  mild  pressure 

_ 1J 4.  / 1 „  4 .U«  _ _ 4-J .. 


gradient  (so  long  as  the  assumption 


-g-«i 

u  Sx 


is  still  tenable)  the  solution  for  7  must  satisfy  the 
non-homogeneous  diffusion  equation! 


g+F(x,y)  =k0- 


with  appropriate  boundary  conditions. 


2.  INITIAL  PROFILES 


The  initial  profiles,  Tq(o,t])  and  u(o,q),  required 
in  equation  (12)  for  evaluating  the  development  of  the 
velocity  and  total  temperature  distributions  were  obtained 
as  follows.  It  was  assumed  that  the  approaching  boundary- 
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layer  (in  the  X,  Y  plane)  satisfies  the  conditions  for 
Blasius  flow.  Hence,  the  velocity  distribution  satisfies: 


where : 


T  is  the  (dimensional)  free-stream  static  temperature. 
The  inital  profiles  used  in  the  present  study  are 
shown  in  Fig.  (48)  for  the  values  Tw  =  2.058*  and  zero 

o 

*This  corresponds  to  an  Eckert  number,  E  =  U  /C  T  (T  -  1), 
of  2.36.  P  ~  w 
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heat  flux  to  the  wall,  these  being  used  as  representative 
conditions.  As  noted  in  Chapter  V,  hot  wire  measurements 
were  unattainable  in  the  separation  corner  region  due  to 
wire  breakage,  and  the  experimental  conditions  for  the 
initial  profile  cannot  be  modeled  exactly.  The  initial 
profile  was  calculated  for  a  Prandtl  number  of  0.72. 

No  allowance  was  made  for  the  distortion  of  the 
profiles  through  the  expansion  process  at  the  corner,  nor 
was  the  presence  of  the  lip  shock  considered.  These 
approximations  are  justified  as  the  present  aim  is  only 
to  study  the  shear  layer  development  in  a  qualitative 
manner.  The  results  should  be  evaluated  in  this  context. 

3.  NUMERICAL  RESULTS  AND  DISCUSSION 

The  initial  velocity  and  temperature  profiles  used 
in  the  study  are  exhibited  in  Fig.  (48).  The  static  tempera¬ 
ture  gradient  at  the  wall  is  zero  since  zero  heat  flux  was 
assumed.  For  a  tunnel  total  temperature  of  560°R,  this 
could  be  representative  of  conditions  after  the  tunnel 
has  been  running  for  some  time. 

As  suggested  by  the  result  shown  in  Fig.  (49),  the 
static  temperature  distribution  changes  very  slowly  along 
the  shear-layer.  It  should  be  noted  that  this  result  is 
for  laminar  flow  and  if  transition  takes  place  in  the  shear 
layer  the  mixing  rate  would  be  somewhat  higher.  The  transi¬ 
tional  case  has  not  been  calculated  due  to  lack  of  data 
concerning  an  acceptable  shear  stress  model  for  such  flows. 

The  same  slow  variation  with  distance  along  the  shear 
layer  is  also  a  property  of  the  velocity  distribution 
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INITIAL  PROFILE 


and  hence  of  the  total  temperature  distribution. 

The  shear  layer  total  temperature  distributions 
are  shown  in  Fig.  50.  The  experimental  data  from  Fig.  46 
is  also  included  on  Fig.  50.  The  comparison  shown  here 
should  be  treated  with  some  caution  since  the  exact 
location  of  the  dividing  streamline  is  not  known  in  the 
experimental  case.  In  Fig.  50,  the  dividing  streamline 
location  indicated  in  Fig.  46  was  used.  Since  it  is  to 
be  expected  that  the  initial  profile  has  its  minimum 
total  temperature  at  the  wall  (i.e.  in  the  absence  of 
extreme  heat  transfer  rates  from  the  wall) ,  and  as  the 
(albeit  simple)  theory  just  presented  indicates  very  small 
rates  of  change  of  temperature  along  the  shear  layer,  it 
can  be  expected  that  the  dividing  streamline  lies  close 
to  the  minimum  total  temperature  condition.  While  it  is 
true  that  the  theory  has  neglected  the  flow  in  the  cavity, 
with  its  associated  heat  transfers,  there  seems  to  be  no 
evidence  suggesting  that  large  temperature  gradients  along 
the  shear  layer  will  exist  for  laminar  flow.  For  this 
reason,  it  is  felt  that  the  comparison  shown  in  Fig.  50 
is  qualitatively  correct. 


Fig.  50.  Comparison  with  experimental  data  (from  Fig.  46) 
(Tentative  only  since  origin  is  unknown  for  the 
experimental  data) . 


CHAPTER  X 


THE  WAKE  SHOCK 

The  separated  flow  region  is  closed  by  the  reattach¬ 
ment  of  the  flow  to  the  downstream  surface.  Since  this  reattach¬ 
ment  process  involves  a  compressive  turning  cf  the  flow  into  the 
free-stream  direction,  a  shock  wave  (the  wake  shock)  is  gen¬ 
erated  at  reattachment.  The  strength  and  inclination  of 
this  shock  wave  is  dependent  upon  its  upstream  conditions 
and  the  extent  of  the  turning  that  is  required.  One 
characteristic  of  a  supersonic  base  flow  is  that  the  length 
of  the  free  shear  layer  decreases  as  the  unit  Reynolds 
number  increases*  Hence,  it  follows  that  the  Mach  number 
along  the  shear  layer  (Mshear)  will  increase  with  increasing 
Reynolds  number  -  and  the  flow  conditions  ahead  of  the  wake 
shock  are  dependent  upon  the  Reynolds  number. 

Experimental  data  [26]  has  indicated  a  very  slight 
sensitivity  of  the  wake  shock  angle  to  changes  in  Reynolds 
number,  which  would  appear  contradictory  to  the  above  remarks. 

The  summarized  experimental  data  [26]  is  shown  in  Fig.  51 

g 

for  a  Reynolds  number  range  of  0.5  -  2.0  x  10  .  - 

In  this  chapter  a  simple  flow  model  is  put  forward 
which  demonstrates  -  in  regions  where  the  viscous  forces 
are  not  significant  -  the  pertinent  features  of  such  base 
flows.  It  is  of  interest  to  know  both  the  reattachment 

*or,  equivalently,  the  free  stream  static  pressure  for 
constant  total  temperature. 
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point  location  and  the  wake  shock  inclination  so  that  the 

/ 

associated  heat  transfer  rates  and  the  far  field  shock 
interaction  phenomena  can  be  designed  for. 

2.  FLOW  MODEL 

The  development  of  the  flow  model  is  the  following. 
Firstly,  the  complex  flow  in  the  vicinity  of  the  corner 
(see  Refs.  42,  48)  is  replaced  by  an  equivalent  simple 
expansion  process.  The  magnitude  of  this  expansion  being 
such  as  to  expand  the  flow  to  the  experimentally  determined 
base  pressure.  In  the  present  work,  the  base  pressure  was 
determined  from  the  empirical  correlation  presented  by 
Su  and  Wu  [46]  and  is  outlined  in  Chapter  VI.  Such  a 
simple  flow-model  omits  the  details  of  the  viscous  transonic 
flow  close  to  the  corner  and  also  ignores  the  presence 
of  the  lip  shock  (which  is,  of  course,  a  product  of  the 
corner  flow) . 

The  comparison  data  shown  in  Fig.  52  indicates,  however, 
that  the  flow  model  predicts  adequately  the  turning  angle 
and  hence  the  length  of  the  shear  layer.  This  is  important 
since  the  peak  heat  transfer  rate  occurs  just  downstream 
of  the  reattachment  point  [25]  and  so  its  location  can.  be 
predicted.  The  flow  properties  ahead  of  the  wake-shock  are 
then  known. 

Conitions  downstream  of  the  wake  shock  can  be  specified 
in  one  of  two  ways;  either  in  terms  of  the  overall  pressure 
ratio,  or  as  a  function  of  the  downstream  flow  direction. 
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In  either  case,  the  downstream  condition  is  not  known,  and 
some  empirical  data  must  be  used.  Herein  it  has  been  as- 
summed  that  the  flow  is  directed  along  the  wall  downstream 
of  the  wake  shock. 


3.  DISCUSSIONS 

The  flow  model  can  now  be  discussed  to  determine  the 
reason  for  the  apparent  invariance  of  the  wake-shock  angle 
with  change  of  Reynolds  number*,  over  the  range  0.5  -  2.0 
x  10  of  the  latter  variable.  Fig.  53  indicates  that 
the  rate  of  change  of  base  pressure  decreases  with  increase 
of  Reynolds  number  over  the  Reynolds  number  range  of 
interest,  and  this  base  pressure  variation  produces  the 
resultant  shear  layer  Mach  number  (MShear)  sJlown  in  53  ♦ 

The  variation  of  the  shear  layer  inclination  was  implied 
in  Fig.  52.  It  is  clear  from  Figs.  52  and  53  that 
the  rate  of  change  of  upstream  wake-shock  variables  is 
still  too  rapid  to  account  for  the  constancy  of  the  wake- 
shock  inclination  (Fig.  51),  and  some  additional  phenomenon 
must  preside. 

It  is  noted  that  the  inclination  of  the  shock-wave  to 
the  free-stream  direction  has  two  components.  One  of  these 
is  the  inclination  of  the  shock-wave  to  the  shock-upstream 
flow  direction  (0),  the  other  is  the  inclination  of  the 
shear  layer  to  the  free  stream  direction  (v).  The  notation 

♦Reynolds  number  (Re)  is  based  on  length  to  the  corner. 
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Fig.  53.  Properties  of  the  flow  model. 


is  shown  in  Fig. 54  .  A  plot  of  6  vs  v  from  the  flow 
model  is  presented  in  Fig.  54,  and  this  clearly  shows 
that  in  the  Reynolds  number  range  of  interest,  the  rates 
of  change  of  both  6  and  v  are  almost  identical.  They  are 
identical  on  the  lines  marked  f3  constant  -  which  the  model 
approaches  at  the  higher  Reynolds  numbers.  Thus,  it  is 
suggested  that  the  constancy  of  the  wake-shock  inclination 
results  from  the  fact  that,  for  a  certain  Reynolds  number 
range,  the  rate  at  which  the  wake  shock  rotates  with  re¬ 
spect  to  its  upstream  direction  is  almost  the  same  as 
that  at  which  the  upstream  flow  is  rotating  with  respect 
to  the  free  stream  direction.  This  contingency  is  the 
result  of  both  the  form  of  the  base  pressure  variation  shown 
in  Fig.  53  and  of  the  oblique  shock-wave  surface  geometry. 

The  results  can  be  extended  to  other  configurations 
(e.g.  wedge  type  flows  discussed  in  Ref.  8)  if  an  appro¬ 
priate  base  pressure  relation  is  introduced. 


CHAPTER  XI 


CONCLUSIONS  AND  RECOMMENDATIONS 

A  detailed  flow  study  has  been  conducted  in  the 
regions  upstream  and  downstream  from  the  separation  corner 
of  a  rearward  facing  step  at  supersonic  speeds.  The 
important  results  and  conclusions  of  this  study  are  pre¬ 
sented  below. 

A  more  direct,  accurate  and  simple  hot-wire  data 
reduction  method  was  developed.  The  convective  heat  trans¬ 
fer  rate  and  the  adiabatic  wire  temperature  were  directly 
determined  from  the  measurements  without  introducing  the 
end  loss  correction  factors.  This  method  makes  possible 
successful  hot-wire  measurements  in  the  free  shear  layer, 
and  in  the  low  Reynolds  number,  subsonic,  recirculating  flow. 

The  experimental  information  concerning  the  upstream 
influence  of  the  lower  base  pressure  is  presented.  A 
significant  fraction  of  the  surface  pressure  drop  occurs 
immediately  upstream  of  the  separation  corner.  The  mag¬ 
nitude  of  the  upstream  influence  is  a  function  of  the  thick¬ 
ness  of  the  undisturbed  boundary  layer.  Correlation  of 
the  measured  surface  pressure  with  the  undisturbed  boundary 
layer  thickness  indicates  that  the  effect  of  the  lower  baso 
pressure  is  propagated  upstream  to  a  distance  of  about 
five  boundary  layer  thicknesses  ahead  of  the  corner  and 
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that  the  sonic  line  is  very  close  to  the  wall  at  the 
separation  corner. 

By  a  simple  order  of  magnitude  analysis,  one  finds 
that  the  acceleration  and  expansion  of  the  major  portion 
of  the  oncoming  boundary  layer  over  the  step  may  be  approx¬ 
imated  as  an  inviscid  process.  An  isentropic  streamtube 
analysis  of  the  measured  pitot  profiles  in  the  approaching 
boundary  layer  shows  that  considerable  acceleration  takes 
place  in  the  lower  half  portion  of  the  approaching  bound¬ 
ary  layer;  and  that  the  normal  and  axial  pressure  gradients 
near  the  corner  are  of  about  the  same  magnitude.  The  static 
pressure  variation  in  the  supersonic  portion  of  the  boundary 
layer  is  roughly  linear. 

The  measured  pitot  profiles  of  the  initial  free 
shear  layer  indicate  that  the  rapid  corner  expansion  cannot 
be  described  by  the  Prandtl-Meyer  theory,  when  the  boundary 
layer  at  separation  is  of  finite  thickness.  It  also 
indicates  that  the  oncoming  boundary  layer  has  a  signifi¬ 
cant  effect  in  the  development  of  the  free  shear  layer. 

The  interaction  of  the  lip  shock  with  the  free  shear  layer  may 
be  more  important  at  higher  free  stream  Mach  number  and  low 
Reynolds  number. 

Comparison  of  measured  pitot  profiles  with  the  numeri¬ 
cal  results  from  the  method  of  rotational  characteristics 
appears  to  indicate  that  the  inviscid  rotational  theory 
is  adequate  for  predicting  the  expansion  of  the  supersonic 
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boundary-layer  at  the  separation  corner;  and  also  the 
lip  shock  formation.  It  was  found  that  the  pitot 
pressure  probe  was  not  sensitive  enough  to  detect  the 
origin  of  the  lip  shock.  The  actual  location  of  the 
lip  shock,  and  its  region  of  formation,  are  a  function 
of  the  Reynolds  number  and  free-stream  Mach  number.  In 
some  instances  the  lip  shock  may  form  some  distance 
downstream  of  the  separation  corner. 

It  was  found  that  the  parameter  Re0*^  (L/H)0*1  was 

oo ,  H 

satisfactory  for  correlating  the  base  pressure  when  the 
shear  layer  flow  was  either  laminar  or  transitional.  The 
limited  amount  of  testing  carried  out  on  the  step  with 
overhang  was  inclusive,  but  appears  to  indicate  that  the 
overhang  has  little  effect.  Further  study  should  be 
carried  out  in  this  area. 

Theoretical  studies  on  the  temperature  distribution 
in  the  shear  layer  indicate  that  the  detailed  temperature 
distribution  in  the  shear  layer  is  most  strongly  influenced 
by  the  initial  profile;  the  effects  of  mixing  in  the  layer 
being  small. 

The  basic  features  of  the  recirculating  region  are 
illustrated  in  terms  of  Mach  number  and  total  temperature 
mappings.  The  dividing  streamline  and  minimum  velocity 
lines  are  identified.  The  total  temperature  of  the 
recirculating  region  may  be  as  low  as  0.8  of  the  free 
stream  total  temperature,  while  the  Mach  number  in  the 
recirculating  flow  may  be  as  high  as  0.8. 


A  simple  model  of  the  inviscid  flow  was  put  forward  to 
study  the  observed  fact  that  the  wake  shock  angle  was  only 
a  weak  function  of  Reynolds  number.  It  was  shown  that  the 
model  was  quantitatively  correct  over  the  Reynolds  number  range 
of  the  test,  but  indicated  a  larger  variation  of  wake  shock 
angle  outside  this  range. 

Based  on  this  study,  the  following  recommendations 
are  presented. 

The  newly  developed  hot-wire  method  is  ideal  for 
measuring  h  and  T  in  the  flow  field  when  the  end  loss 
corrections  are  required.  This  technique  is  recommended 
to  check  the  existing  heat  transfer  correlation  formula 
and  if  necessary,  to  provide  new  heat  transfer  data. 

A  heat  transfer  correlation  formula  is  needed  for 
low  Reynolds  number  and  very  low  subsonic  flow  (ReQ  <  0.20, 

M  <  0.3).  Experimental  information  concerning  natural  convec¬ 
tion  in  this  range  is  also  needed.. 

It  is  desirable  to  verify  the  estimation  of  turbu¬ 
lent  effects  on  measured  hot-wire  resistance  and  to  meas¬ 
ure  accurately  the  turbulent  intensity  in  the  free  shear 
layer  and  the  recirculating  flow.  Hot-wire  measurements 
in  the  recirculating  region  with  laminar  free  shear  layer 
are  needed  for  comparison  with  theoretical  results. 

Many  free  shear  layer  profiles  are  presented  at 
different  Mach  numbers,  Reynolds  numbers  and  step  heights. 
Sufficient  experimental  information  is  provided  for  theo¬ 
retically  analyzing  the  problems  related  to  the  free  shear 
layer  growth. 
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Both  theoretical  analysis  and  accurate  and  detailed 
measurements  are  needed  to  study  the  acceleration  and 
expansion  of  the  subsonic  portion  of  the  oncoming  boundary 
layer,  which  hopefully  will  completely  define  the  initial 
conditions  of  the  mixing  region. 

At  higher  free  stream  Mach  number,  the  shear  layer 
is  expanding  and  mixing  with  the  recirculating  flow;  and 
the  effects  of  pressure  gradients  on  the  mixing  process 
deserve  attention. 

A  comparison  of  the  hot-wire  results  with  numerical 
calculations  of  mixed  supersonic/subsonic  type  flow  in 
the  near  wake  [17]  may  be  helpful  for  interpretation  of 
the  hot-wire  measurements.  The  successful  numerical  solu¬ 
tion  of  the  compressible  Navier-Stokes  equation  is  very 
difficult,  however,  successful  measurements  in  part  of  the 
near  wake  (such  as  separation  region  and  re attachment 
region)  are  also  very  difficult.  One  should  not  rule  out 
the  potential  of  the  numerical  method  in  solving  the  near 
wake. 

The  upstream  influence  was  studied  by  Weiss  and 
Nelson  [5]  with  no  static  pressure  gradient  in  the  accel¬ 
erating  boundary  layer.  It  may  be  interesting  to  follow 
their  approach  but  with  linear  static  pressure  distribu¬ 
tion  in  the  supersonic  portion  of  the  boundary  layer. 
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APPENDIX  A 


REVIEW  ON  HOT-WIRE  MEASUREMENTS  IN  LOW  REYNOLDS 
NUMBER  HYPERSONIC  FLOWS 


The  technique  of  hot-wire  measurements  in  low 
Reynolds  number  hypersonic  flows  was  originally  dev<- Toped 
by  Dewey  [31].  The  technique  is  described  as  follows. 

The  convective  heat  transfer  from  length  Z  of  a 
heated  wire  of  infinite  length  to  the  surrounding  fluid 
is 


lj|  R*  =  hrrdfc  (T*  -  T^)  (A-l) 

where  superscript  *  refers  to  infinitely  long  wire  and  all 
the  assumptions  made  in  Chapter  III,  Section  2  are  re¬ 
tained.  The  term  "infinite  wire"  refers  ;o  a  wire  of 
sufficient  length  that  end  loss  is  not  important.  It 

should  be  emphasized  that  both  h  and  T  are  assumed  to 

aw 

be  independent  of  the  wire  length. 

From  Equation  (A-l),  the  convective  heat  transfer 
coefficient  is 


h  = 


Zm  R* 
m 


wdMT*  -  T  ) 
aw 


(A- 2) 


For  a  finite  wire,  the  measured  convective  heat 
transfer  coefficient  is  defined  as 
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11  R„  =  h„Trd£(T  -  T,w  ) 


mm  m 


hm  “ 


I2  R 
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*dll<Tm  '  Taw  ) 
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(A- 3) 


where  the  subscript  m  denotes  quantities  measured  by  a 

finite  wire,  and  the  heat  transfer  includes  both  convection 

and  conduction.  In  order  to  correct  h  to  the  value  h 

m 

for  an  infinitely  long  wire,  one  introduces  the  ratio  of 
these  two  quantities  at  a  given  current.  This  ratio  is 
defined  as  the  Nusselt  number  correction  factor 


h  Nu  R  T  ~  Taw 

^  _  11  _  _ o  _  ___  m  awm 

N  h  Nu  R  T*  -  T_w 

m  m  m  aw 


(A- 4) 


Similarly  the  correction  factor  for  the  recovery 


factor  is  defined  as 


h  = 


(A-5) 


Tawm  ^ 


i^N  and  \pR  are  derived  from  the  measurements  with  a  finite 
wire  in  the  following  way. 

From  Equation  (11),  the  measured  wire  temperature 


1/2 

2  f  n  [(n/m2)  -  T  ] 

T  =  —  I  T  dx  =  — j  -  -  tanh 

m  Z  J  m*  mi/2 


(A- 6) 
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where 


mil  hd  b  I2  °*5  i 

_  /  m. 


2  kw  d 


for  an  infinitely  long  wire, 


i 

d 


mil 

(tanh  — )/{mil/2)  0 

2 


and 


* 


T 


lim 
%/d  * 


n 


(A- 7) 


bn 


R  =  a  + 


(A- 8) 


m 


Substituting  Equations  (A-7)  and  (A-8)  into  Equation 
(A-4),  one  obtains 


a  +  (bn/m  )  T  -  Taw 

=  £— ; - 1  t- St  —  -) 


(A- 9) 


m 


(n/m  )  -  T 


aw 


After  a  lengthy  expansion  and  neglecting  the  higher 
order  terms  with  some  approximations,  Dewey  finally 
obtained 


%  -  V*>  ■  *H(NV 


(A-10) 


since 


l  hd  0-5  l  k  °*5 

S  =  —  ( — )  =  -  (Nu  — ) 

d  k  d  °  k 

w  w 

Substituting  Equation  (A- 7)  into  Equation  (A- 6) , 
one  obtains 

tanh  (mi/2)  tanh  (m£/2) 

Tm  =  T*  [1 - ]  +  Ts[ - ]  (A-ll) 

■(mJl/2)  (m£/2) 

and 


Tawm  " 


T  =  T  {1  - 
m  aw 


tanh  S  tanh  S 

- )  +  t  ( - ) 


I  0 


(A-12) 


where 


*aw  T 

I  -*■  0 


Rearranging  Equation  (A-12) ,  one  gets 


or 


*R  = 


aw 


T 


aw 


1  -  (Ts/Tawm)  (tanh  S/S) 
[1  -  (tanh  S/S)] 


m 


(A-13) 


’J'r  =  «Pp(S)  =  ’I'd  (Nu„) 


R 


(A-14 ) 
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Then,  Equation  (A-10)  and  (A-14),  and  the  correlation 
formula  for  heat  transfer  and  recovery  factor  were  coupled 
together  to  determine  the  flow  field  parameters,  i.e., 

Taw(M,  xn>  =  TaWm 

NU0(M,Re0)  =  Num  *R(Reo) 

The  method  of  successive  substitution  was  used  to 
solve  the  above  system  of  equations  and  convergence  of 
the  solution  was  not  obtained  when  M  <  1.5. 

At  higher  Mach  number  (M  >  2),  the  correlation 
formula  for  Taw  and  h  becomes  more  simple.  Behrens  [52] 
used  the  following  two  equations  for  the  data  reduction  of 
his  hot-wire  measurements 

Ite  »  5.26  Nu_  +5.74  Nu2 
o  o  o 


T 


o 


1.167 


0.217 

1  +  [1/(0.335  Reo)1,33] 


136 


APPENDIX  B 


SOLUTION  FOR  h  AND  T 


aw 


Equations  (13)  and  (14)  are  solved  for  h  and  T 


aw 


by  a  Newton-Raphson  iteration  scheme.2  The  details  are 
described  as  follows.  Equations  (13)  and  (14)  can  be 
rearranged  as 


b(Taw  "  Tjtanh  S 

F(h,Taw)  =  (a  +  b  Taw)  -  — SS. - f - Rm 


(B-l) 


I  -*•  0 


a  +  b  T 


G(h,T  )  =  b(- 


aw 


tanh  S 


aw 


irdWi 


)  (1  - 


-) 


b2  (T  -  T  )  sech2  S  tanh  S  dR 

+  — 3S - SI  ( - - - )  - 


2S 


irdtdidk  ,  irdJlh  dl 
w  m 


(B-2) 


|I  0 

The  partial  derivatives  of  F  and  G  with  respect  to 


h  and  T  .  are 
aw 


3F 


3h 


S  sech  S  -  tanh  S 


"  -  b(Taw  -  Ts)  (- 


2h  S 


-) 


(B-3) 


3F  tanh  S 

- —  =  b  (1  - - ) 


( B—  4 ) 


3T 


aw 


^The  author  is  indebted  to  Dr.  K.  C.  Reddy  for 
suggesting  this  method. 


a  +  b  T 


tanh  S  S  sech  S  -  tanh  S 


-)  (  1  - 


b2(T  -  T  )  sech2  S  tanh  S 
_ aw _ s _  ^ _  _  _ 

4hSird  ✓hdJc  £h 

w 


b2(T  "  Ta>  S  (sech  S)2  tanh  S 

QW  S  ^ 


2hSnd 


^aig- 


2 

sech  S 


2^ 


S  sech  S  tanh  S 

+ - ] 

2£h  £h 
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3G  b 


3T  ird£h 

aw 


2 '  2 

tanh  S  b  sech  S  tanh  S 

(1 - )  +  __  ( - ) 

S  2S  irdvfiaF"  Trd£h 


(B-6) 


Equations  (B— X)  through  (B-6)  are  then  substituted  into 
the  following  two  equations,  which  are  the  Taylor  series 
approximations  of  F  and  G. 


F(VTaWi)  =  Ahi  “  +  ATaw.  “ 

1  i  1  3h  .  i  3T 


aw  . 

l 


(B-7) 


G(hi,Taw.)  =  Ahi  --  +  ATaw_  — 

i  3h  .  l  3Taw  i 


(B-8) 


Solving  the  above  two  equations  for  Ah^  and  ATaw  ,  then 

th  i 

the  improved  values  for  (i  +  1)  iteration  are 
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APPENDIX  C 


SOLUTION  FOR  M  AND  T 

o 


The  Newton-Raphson  iteration  scheme  is  also  employed 

for  the  numerical  calculation  of  M  and  T  .  The  details 

o 

are  described  as  follows.  Equations  (16)  and  (17)  can 
be  rearranged  as 


F(TQ,M) 


hd 

- 5(Re0)  iMRe0,M) 

k 


(C-l) 


T 

G(To,M)  =  —  -  nc  -  n*  (nf  -  nc) 
T  o 


(C-2) 


Then  the  approximated  equations  for  F  and  G  at  the  i 
iteration  are: 

(C-3) 

(C-4) 

The  partial  derivatives  of  F  and  G  with  respect 

to  T  and  M  are: 
o 


8F 

F(Tq  ,M.  )  =  ATq  - 

i  1  i  3T. 


+  AM. 


3F 


3M 


3G 

G  (Tq  ,M.)  =  ATq  — 
i  1  i  3T, 


+  AM. 


3G 


3M 
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3F 

3T 


0.86  hd 


3? 


dip  3Re 


6.7  x  io"5  x  T*'86 
o  o 


-  (t - 5 - ) 

3Re  3Re^  3T 

o  o  o 


(k  =  6.7  x  io"5  T*86  BTU  ft-1  nr-1  0R-1) 
o  o 


3F  3?  3Re  dip 

- -  ♦ - -  -  S  — 

3M  3Re  3M  3M 

o 


3G  Taw  0.589 (nf  “  0C)  Kn°‘193  SKn^ 


3T 


T 


(0.493  +  Kn**193)2  3Tq 


3G 

3M 


.2.5 


0.206  o.589(nf  -  nc) 


Kn 
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3Kn 


(1.175  +  M3*5)2 


(0.493  +  Kn*’193)2  3M 


0.516  rT*  M1*8 
(0.8521  +  M2,8)2 
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H 

3Re. 


=  (0.14  +  0.230  (- 


Re 


0.7114 


15.44  +  Re 


0.7114 


) 


+  ( 


0.01569 


F>  (■ 


15 


0.3077  +  Re0,7378  15  +  Re3 

o  o 


)] 
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nrl  .  0.862 (In  Rep)  Re^ 

‘  Reo  Re0.3287  (2  571  +  ^Q.6713,21 


n  2.524  Re"0,2886 

+  Re»^  f  ^ 

°  (15.44  +  Re°‘7114)2 


0.7065  Re* 

(0.3077  +  Re8,7378)  (15  +  Re2)2 


0.1737  Re“°*2622 

(15  +  Re3)  (0.3077  +  Re0,7378)2 
o  o 


5.01  Re9,109  0.065  Ren 

- Vl-Qo-y  U  +  (0.30  -  -™-)  ( - 2__)] 

(2.765  +  Re^  uyr  Mx'b/  4  +  Re 
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+  [1.834  -  1.634  (- 


2.675  +  Re' 


. .  109 
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x  [(0.30  -  ( - 

M1*67  4  +  Re 


2  2 
- )  ] 


0.6039  1.0929  M 

- —  +  - - — — 
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3Re, 


M2  (1  +  m1*222)2’569  9M 


Re?;'109  Re0  0.1085 

+  [1.834  -  1.634  ( - 2 - __)  ]  [( - )  (“TW 

2.675  +  Re1,109  4  +  Re  M2,67 


APPENDIX  D 


STATIC  PRESSURE  DISTRIBUTION  DOWNSTREAM  OF  BASE 

The  surface  pressure  distributions  downstream  of 
the  base  are  presented  in  Figures  55 ,  56  and  57  for  three 
typical  Reynolds  numbers  at  each  of  the  three  Mach  numbers 
M^  -  2.5,  3.5  and  5.07.  p  is  the  measured  static  pressure, 
p^  is  the  static  pressure  measured  2.5  inches  from  the 
leading  edge,  and  x  is  the  distance  downstream  from  the 
base.  The  boundary  layer  before  separation  was  always 
laminar  except  at  the  conditions  M^  =  2.5,  Re^  L  =  1.96 
x  1C  .  The  flow  field  was  entirely  laminar  until  after 
reattachment  at  the  lowest  Reynolds  number.  At  the  largest 
and  intermediate  Reynolds  numbers,  the  transition  took  place 
somewhere  in  the  free  shear  layer,  and  the  reattachment 
was  therefore  turbulent. 

In  the  pressure  distribution,  tnere  is  a  hump  in 
the  midst  of  pressure  recovery  process  and  a  local  peak 
after  the  pressure  is  sufficiently  recovered.  This  hump 
and  local  peak  are  attributed  to  the  interaction  of  the 
lip  shock  with  the  wake  shock  [8,  20].  However,  it  is 
found  that  this  hump  and  peak  may  exist  simultaneously 
under  certain  conditions. 
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Figure  55 .  Pressure-Recovery  Distribution 


Figure  57.  Pressure-Recovery  Distribution  at  M, 


APPENDIX  E 


EFFECT  OF  VELOCITY  FLUCTUATION  ON 
RESISTANCE  MEASUREMENTS 

To  simplify  the  analysis,  the  end  losses  due  to  heat 
conduction  are  neglected.  The  heat  balance  equation  is 
then 


ll  R  =  hA(T  -  T  )  (E-l) 

mm  m  aw 

The  convective  heat  transfer  formula  is  also  assumed  to 
take  the  simplified  form 

Nuo  -  °1  +  °2  ReS  (E-2> 

Since  only  the  velocity  fluctuation  is  considered. 
Equation  (E-2)  becomes  [53] 

h  -  c3  +  c4  «e£f  <E'3> 

and 

Ueff  =  &  +  u')2  +  v'2  +  w'2  (E_4) 

where  U  is  the  steady  velocity  and  u',  w' ,  w",  the  fluc- 
uation  velocity  component  orthogonal  to  one  another. 
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The  instantaneous  wire  resistance  is 

R  =  R  +  r"  (E-5) 

m  m 

Where  R  is  the  steady  resistance  and  r'  is  the  fluctua¬ 
nt 

tion  resistance  due  to  velocity  fluctuation. 

The  linear  relationship  between  temperature  and 
resistance  is  assumed 

R  *  a  +  b  T  (E-6) 


Substituting  Equations  (E-3)  and  (E-6)  into  Equation  (E-l) , 
one  obtaines 


or 


■  ~  <*»  '  Raw>  <°5  +  °6  Ue£f> 


(E-7) 


I2  (R_  +  r")  »  -  (R  +  r' 
mm  b  m 


-  *,J  (=5  +  c6  U”£f) 


aw' 


(E-8) 


By  series  expansion 


2u'  m'2  +  v'2  +  w'2  n/2 

U 


Ueff  =  [U2  (1  +  —  +  - — - )] 


U2 


=  un  (1  +  — -  +  -  (- 
U  2 


nu'  n  u"2  +  v"*2  +  W*2 


u2 


•)  +  ...  3 

(E-9) 


Substituting  Equation  (E-9)  into  (E-8),  taking  the  mean 
time  average  and  neglecting  the  higher  order  terms. 
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one  obtains 


_n  n  2 

Cg  fr  (1  +  -  T^)  +  C5 
6 


aw  c.  Un  (1  +  -  T2)  +  c.  -  I2  b 
6  6  5  in 


(E-10) 


where 


T  m 


[-  (u'2  +  v'2  +  w'2)J 
3 


0.5 


U 


which  is  the  intensity  of  turbulence  of  a  stream. 

Expanding  Equation  (E-10)  and  neglecting  the 
higher  order  terras,  one  gets 


R 


ra 


aw 


( 


c6  u”  +  c5 


c,  Un  +  cK  -  I2  b 
6  5  m 


n 

HI-- 


T2  cc  Un  I2  b 
o  ra 


6  <c5  Un+o5)(c6  Un+  c5-  IjJ  b) 


] 


(E-ll) 


Bauer  obtained  the  same  expression  using  a  different 
approach  [54 ] . 

From  Equation  (E-7),  one  knows 


Im  b 


cg  U  +  Cg 


Rra  ~  Raw 
Rra 


also 
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1 


=6^ 


°6  un  +  c5  '  Tm  b 


0.5  <  n  <  1 


If  one  assumes  a  10  percent  velocity  fluctuation 
in  the  shear  layer  and  recirculating  region,  then  the 
percentage  decrease  in  measured  resistance  due  to  turbu¬ 
lence  is 


n  c.  Un 

(-)  ( —zr1 - *->  <- 


e  o,  u"  +  c5  -  i;  b  c6  s"  +  c5  -  i;  b 


I2  b  2 

— - o— )  T2 


Raw  2 

=  -  (1  -  — )  T  s 


(.04)  (.01) 


=  .000067  «  1 


6  R 


m 


Therefore,  one  may  conclude  that  the  effects  of  turbulence 
on  the  resistance  measurement  is  negligibly  small. 
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APPENDIX  F 


REMARKS  ON  HOT-WIRE  MEASUREMENTS 

1.  HOT-WIRE  DIAMETER  AND  UNCERTAINTIES 

As  discussed  in  Chapter  III,  it  is  very  difficult 
to  apply  the  conventional  calibration  technique  for  hot¬ 
wire  measurements  in  low  Reynolds  number  compressible 
flow.  However,  the  technique  which  relies  completely 
upon  the  empirical  heat  transfer  correlation  formula  is 
by  no  means  perfect.  There  are  inherited  difficulties 
and  uncertainties  associated  with  this  technique. 

The  diameter  of  the  hot-wire  is  about  0.0002  inch. 

It  cannot  be  made  exactly  circular  and  uniform.  The  cross- 
section  of  the  wire  may  be  elliptical  or  irregular,  and 
there  is  no  unique  way  of  determining  the  wire  diameter. 

If  one  assumes  that  the  wire  is  circular  and  of  uniform 
cross-section,  then  the  wire  diameter  can  be  determined 
by 

pi 

R  =  — 2— 

TTdV4 

where  R  is  the  total  resistance  of  the  wire  of  length  £ 

2 

and  cross-sectional  area  ird  /4,  and  p  is  the  resistivity 
of  the  materiel. 
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The  wire  diameter  was  estimated  to  be  0.000217  inch 
by  the  above  equation  and  with  the  following  information 
(p  supplied  by  wire  manufactueres) . 

I  =  0.1354  ±0.0002  inch 
R  =  26.57  ohm 
p  =  7.261  *  10~®  ohm  inch 

Since  the  hot-wire  diameter  is  about  seven  wave¬ 
lengths  of  visible  light,  optical  methods  for  determining 
the  wire  diameter  are  not  very  accurate. 

The  physical  constants  of  the  materiel  are  not 
available  for  extremely  thin  wires.  For  instance,  the  wire 
thermal  conductivity  may  be  different  from  that  of  the 
same  materiel  with  larger  dimensions.  This  may  introduce 
some  error  in  calculating  the  conduction  heat  transfer 
to  the  hot-wire  supports. 

The  reduced  results  of  the  hot-wire  measurements 
depend  entirely  on  the  accuracy  of  the  correlation  formulas. 
These  correlation  formulas  are  not  likely  to  be  suffi¬ 
ciently  accurate  at  the  hypersonic  and  transonic  Mach 
numbers.  It  is  also  very  difficult  to  obtain  satisfactory 
correlation  formula  at  low  Reynolds  number  and  in  the  low 
subsonic  range,  since  the  data  in  this  regime  are  very 
limited  and  diverged. 
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2.  RADIATION  HEAT  LOSS 


The  radiation  heat  transfer  from  the  hot-wire  was 
neglected  in  the  data  reduction/  and  it  remains  to  be 
justified.  The  ratio  of  the  radiation  heat  loss  to  the 
convective  heat  transfer  is 


a(T* 
_ m 

h(Tm 


) 


T 


aw 


) 


where  a  is  the  Stef an- Boltzmann  constant. 

The  calculated  results  given  below  show  that  even 
at  the  low  Nusselt  number  radiation  from  the  hot-wire  is 
negligible. 


M 

Re^ 

o 

Nu 

O 

h 

V#R 

Taw'°R 

V^h 

3.4 

4.4 

0.544 

468 

590 

582 

0.0073 

1.2 

1.06 

0.214 

172 

580 

553 

0.0088 

0.34 

0.27 

0.155 

125 

570 

521 

0.0017 

where  h  is  in  BTU/hr  ft2  °F  and  Tq  =  560°R. 

3.  NATURAL  CONVECTION 


Along  the  minimum  velocity  line  in  the  recircu¬ 
lating  region/  the  velocity  may  be  very  low  and  the  velocity 
induced  by  buoyancy  may  be  of  equal  magnitude.  The  natural 
convection  shall  be  considered  in  the  heat  transfer  analy¬ 
sis.  However,  as  indicated  in  Chapter  III,  when  the  Mach  num¬ 
ber  was  less  than  0.3,  no  converged  solution  was  obtained. 
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When  the  Mach  number  is  greater  than  0.3,  the  velocity  induced 
by  buoyancy  is  small  compared  to  the  total  velocity.  There¬ 
fore,  it  is  reasonable  to  neglect  the  gravity  effect  in 
the  analysis  of  forced  convection  from  the  heated  wire. 

4.  COMPARISON  OF  MEASURED  RESISTANCES 

Figure  51  shows  the  resistance  traces  for  three 
different  currents  at  a  typical  station.  The  comparison 
can  be  better  explained  by  the  following  equation  for  the 
measured  resistance,  i.e., 

a  +  b  (nTQ/^R) 

Rm  ’  1  ‘  <*n  *2  b'h  A> 

which  is  derived  from  the  hot-wire  heat  transfer  equation 

with  the  conventional  end  loss  correction  factors. 

The  above  equation  shows  that  the  effect  of  h  on 

2 

the  measured  resistance  is  amplified  by  Im.  At  very  small 
current,  the  convective  heat  transfer  has  practically  no 
influence  on  the  wire  resistance;  and  the  wire  measures 

essentially  the  recovery  wire  temperature.  At  the  higher 

.  .  2 
heating  current,  the  relative  contributions  of  h  and  I 

on  Rm  depend  on  their  magnitudes.  Note,  the  changes  in 

resistance  traces  when  the  current  is  increased  from  1  to 

7  and  to  9  milliamps  as  indicated  in  Fig.  58. 
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5.  RECOVERY  FACTOR  IN  TRANSITIONAL  REGIME 


The  correlation  of  recovery  temperature  was  dis¬ 
cussed  in  detail  by  Dewey  134],  The  recovery  temperature 
was  correlated  by  Equation  (17)  in  terms  of  normalized 
recovery  ratio  n*/  Mach  number  and  Knudsen  number.  The 
recovery  factor  in  the  transitional  regime  was  computed 
based  on  Equation  (17),  and  is  shown  in  Figure  59. 
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APPENDIX  G 


ORIFICE  DIMENSIONS 


The  detailed  orifice  dimensions  in  the  wind  tunnel 
model  as  shown  in  Figure  2,  page  8,  are  presented  in 
Table  II. 
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ORIFICE  DIMENSIONS  IN  WIND  TUNNEL  MODEL 
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